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CHAPTER  1 


INTRODUCTION 


Transmission  lines  are  usually  thought  of  as  being  the  wire  net¬ 
work  connecting  a  system  of  electronic  devices.  These  wires  are  in¬ 
tended  to  carry  signals  (often  digital  signals  at  bit  rates  in  the 
low  megahertz  frequency  range)  from  one  electronic  device  to  another. 
Ideally  the  signal  received  at  the  second  device  should  be  the  same 
as  the  signal  launched  by  the  first  device.  However,  very  often  the 
connecting  transmission  lines  are  immersed  in  an  electromagnetic 
field  (FIGURE  1-1).  This  field  can  be  generated  by  nearby  radiating 
antennas  (such  as  onboard  aircraft) ,  by  other  transmission  lines 
carrying  high  frequency  signals  (this  is  referred  to  as  wire-to-wire 
coupling  or  crosstalk),  or  just  the  high  frequency  energy  that  con¬ 
stantly  bombards  our  environment  (radio,  television,  microwave  com¬ 
munications,  radar).  Both  electromagnetic  theory  and  experience 
tell  us  that  this  field  can  add  signals  to  those  already  being 
carried  by  the  transmission  lines. 

It  is  important  to  know  the  level  of  the  unwanted  signals  (noise) 
on  the  transmission  lines  from  the  standpoint  of  device  design  and 
from  the  system  analysis  standpoint.  That  is,  the  designer  of  an  in¬ 
dividual  electronic  device  must  know  what  part  of  the  incoming  signal 
is  informational  and  what  part  is  noise  in  order  to  set  noise  level 
tolerances.  The  system  designer  must  know  what  happens  to  a  signal 
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between  the  time  it  leaves  one  device  as  an  output  and  is  received 
at  another  device  as  an  input.  Even  though  a  device  may  work  well  in 
an  isolated  test  environment,  it  must  also  function  properly  in  close 
proximity  to  other  devices  in  the  system.  In  short,  given  a  trans¬ 
mission  line  network  and  given  an  incident  electromagnetic  field,  it 
is  important  to  be  able  to  predict  when  and  by  how  much  the  field  will 
affect  the  signals  being  carried  on  the  transmission  lines. 

From  an  electromagnetic  point  of  view,  the  problem  here  is  that 
the  transmission  lines  are  acting  like  receiving  antennas.  The  sur¬ 
rounding  electromagnetic  field  induces  currents  on  the  wires  of  the 
transmission  network.  Unfortunately,  the  accurate  prediction  of  these 
induced  currents  can  be  a  somewhat  difficult  and  expensive  task.  The 
problem  is  real  and  a  computationally  efficient  model  would  have 
practical  applications.  One  such  model  is  the  "transmission  line 
model  with  distributed  sources",  described  in  Chapter  2.  This  model 
is  economical  and  probably  the  most  efficient  method  of  prediction 
known  but  has  not  been  fully  validated. 

The  transmission  line  model  is  based  on  several  key  assumptions. 
The  model  is  devised  with  a  particular  physical  structure  in  mind. 

In  general,  the  structure  is  composed  of  parallel,  thin  wires  with 
wire  length,  separation  between  wires,  and  wire  radius  related  in  the 
following  way.  The  wire  radius  is  assumed  to  be  much  smaller  than  the 
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separation  between  the  wires,  and  the  separation  is  assumed  to  be  much 
less  than  the  overall  length  of  the  line.  The  use  of  the  term  "thin 
wire"  implies  that,  at  relevant  frequencies,  the  radius  is  small 
enough  such  that  currents  flow  only  in  a  longitudinal  path  (along 
the  length  of  the  wire)  and  none  flow  circumferentially  (around  the 
wire).  (See  FIGURE  1-2.)  The  separation  between  wires  is  assumed  to 
be  small  enough  so  that,  at  relevant  frequencies,  it  is  not  a  signifi¬ 
cant  portion  of  a  wavelength.  This  assures  that  the  TEM  mode  will  be 
the  primary  mode  of  propagation  and  that  higher  order  modes  will  be 
cutoff. 

The  transmission  line  is  terminated  with  some  type  of  load  or 
source.  While  the  source  can  be  an  alternating  current  source  and  the 
load  can  be  a  complex  impedance,  either  is  considered  to  be  a  lumped 
parameter.  That  is,  the  termination  is  considered  to  be  "electrically" 
small,  and  the  source  or  load  occupies  no  physical  space  but  is  only 
a  mathematical  parameter.  Take  for  example  one  end  of  a  pair  of 
transmission  lines  that  are  connected  to  each  other  by  a  resistor. 

(See  FIGURE  1-3.)  The  current  leaving  the  top  wire  is  exactly  the 
same  as  the  current  entering  the  bottom  wire.  The  resistor,  value 
R,  serves  only  to  relate  the  voltage  at  the  termination  to  the  current 
entering  the  termination.  The  resistor  is  "electromagnetically  in¬ 
visible"  and  the  termination  is  not  a  part  of  the  structure.  What 
this  means  is  that  currents  induced  by  some  incident  electromagnetic 
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field  will  only  be  induced  on  the  transmission  wires  themselves,  none 
will  be  induced  in  the  load.  Most  of  these  assumptions  are  closely 
related  to  one  another  and  together  give  transmission  line  theory  its 
economy  and  computational  efficiency,  but  they  also  place  restrictions 
on  the  application  of  the  model. 

Another  technique  used  in  solving  for  induced  currents  is  re¬ 
ferred  to  as  the  "method  of  moments",  which  is  discussed  in  Chapter  3. 
This  method  is  drastically  different  from  the  transmission  line  model 
and  was  designed  specifically  for  use  on  a  digital  computer.  The 
method  of  moments  makes  very  few  assumptions  and  is  a  much  more  rigor¬ 
ous  but  expensive  prediction  method  than  that  of  transmission  line 
theory.  With  the  method  of  moments,  the  entire  structure  is  modeled, 
including  the  termination  load.  This  method  solves  for  currents  on 
the  entire  structure,  rather  than  just  at  the  terminations,  by 
enforcing  electromagnetic  boundary  conditions  throughout  the  entire 
structure.  The  only  limit  to  the  accuracy  of  this  method  is  computer 
core  and  computer  time.  In  lieu  of  experimental  data,  two  method  of 
moments  computer  codes  have  been  selected  to  generate  baseline  data 
with  which  the  transmission  line  model  predictions  can  be  compared. 
These  codes  will  be  discussed  further  in  Chapter  3. 

The  transmission  line  method  solves  for  the  terminal  currents  and 
voltages.  This  makes  it  very  well  suited  for  a  practical  application 
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to  the  "real  world"  problems  discussed  at  the  beginning  of  this  chap¬ 
ter.  The  main  objective  here  is  to  predict  how  currents  induced  on 
transmission  lines  by  an  incident  electromagnetic  field  are  coupled 
into  an  electronic  device  through  its  inputs.  So  while  knowledge  of 
the  actual  currents  at  various  positions  on  the  wires  may  be  impor¬ 
tant  in  some  cases,  such  as  antenna  design,  here  we  are  more  concerned 
with  prediction  of  currents  in  the  terminations  produced  by  fields 
incident  on  the  transmission  lines.  The  method  of  moments,  while 
very  accurate,  is  very  costly.  Terminal  current  predictions  at  200 
frequencies  were  calculated  using  the  transmission  line  model  formu¬ 
lation  on  an  IBM  370/165  digital  computer  in  approximately  2  seconds 
of  CPU  time.  Solution  by  the  method  of  moments  using  the  same  compu¬ 
ter  for  only  10  frequencies  took  an  average  of  80  seconds  of  CPU  time. 
This  means  that  the  transmission  line  model  can  produce  solutions  for 
800  frequencies  in  the  same  amount  of  computer  time  required  for  solu¬ 
tion  at  one  frequency  by  the  method  of  moments. 

It  is  obvious  that  the  method  of  moments  is  not  a  practical  solu¬ 
tion  technique  for  a  large  number  of  frequencies.  However,  its  ac¬ 
curacy  makes  it  an  excellent  source  for  baseline  data  to  validate 
the  transmission  line  results.  In  general  the  method  of  moments 
will  break  down  numerically  (a  digital  computer  cannot  provide  numeri¬ 
cally  stable  solutions)  at  low  frequencies.  The  theory  behind  the 
transmission  line  model  breaks  down  at  high  frequencies.  However, 
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both  techniques  should  provide  good  results  over  a  fairly  large  com¬ 
mon  frequency  range.  In  particular,  the  upper  frequency  limitation 
of  the  transmission  line  model  will  be  investigated.  Here,  the 
method  of  moments  solution  should  be  valid  and  accurate. 

In  Chapter  2  the  transmission  line  model  is  described  in  detail 
and  a  solution  is  derived  for  the  terminal  currents.  In  Chapter  3 
the  basic  principles  behind  the  method  of  moments  technique  are  dis¬ 
cussed.  As  mentioned  earlier,  two  computer  codes  have  been  selected 
to  provide  the  method  of  moments  solution.  Although  both  codes 
use  the  method  of  moments  technique,  each  takes  a  slightly  different 
approach  to  the  solution  and  some  of  the  differences  between  the  two 
codes  are  discussed.  The  structure  used  for  the  majority  of  the  data 
generation  and  the  reasons  for  its  selection  are  discussed  in  Chapter 
4.  A  great  deal  of  data  was  collected  for  this  structure,  for  several 
different  cases.  The  validity  of  the  transmission  line  results  is 
explored  in  detail  here.  Several  special  cases  are  considered  in 
Chapter  5.  These  cases  involve  structures  slightly  different  from 
the  structure  considered  in  Chapter  4.  These  modified  structures  help 
to  demonstrate  several  points  better  than  the  basic  structure.  A  very 
important  topic  is  addressed  in  Chapter  6:  the  relationship  between 
the  common  and  differential  mode  currents  and  their  relationship  to 
the  total  induced  current.  The  findings  presented  in  Chapters  4  and 
5  show  an  amazing  accuracy  for  the  transmission  line  model  predictions 
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of  differential  mode  currents  over  a  considerable  frequency  range. 
These  results  demonstrate  the  integrity  of  this  efficient  model  and 
suggest  that  it  can  provide  accurate  predictions  when  applied  pro¬ 
perly  to  cases  which  are  within  its  limits  of  application. 
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CHAPTER  2 


THE  TRANSMISSION  LINE  MODEL 

The  purpose  of  this  chapter  is  to  describe  the  transmission 
line  model  mathematically.  But  first,  the  parameters  which  physically 
describe  the  structure  must  be  defined.  In  Chapter  4  a  two  wire 
transmission  line  is  analyzed  in  some  detail;  therefore,  the  model 
described  here  will  be  specialized  to  two  wires.  (Refer  to  FIGURE 
2-1 (a).)  As  discussed  in  Chapter  1,  the  structure  consists  of  two 

parallel  wires  terminated  at  each  end  with  complex  loads  Z  and  Z  . 

U  L 

The  line  has  length  L,  which  is  defined  as  being  the  distance  between 
terminal  planes.  The  load  is  not  part  of  the  physical  line  and  makes 
up  no  part  of  the  length.  The  distance  between  the  centers  of  the 
wires  is  the  wire  separation,  d.  The  radius  of  the  wires,  r,  has  no 
effect  on  the  measurement  of  the  separation  between  the  wires  because 
of  the  "thin  wire"  assumption,  which  implies  that  the  radius  (diameter) 
is  insignificant  when  discussing  the  location  of  the  wires.  The 
radius  is  only  used  in  the  calculation  of  per-unit-length  capacitance 
and  inductance,  and  the  characteristic  resistance  of  the  line,  which 
are  discussed  later. 

The  orientation  of  the  transmission  line  with  respect  to  the 
Cartesian  coordinate  system  is  shown  in  FIGURE  2-l(b).  The  structure 
lies  in  the  x-y  plane  (z=0) .  The  bottom  wire  starts  at  the  origin 
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and  extends  along  the  positive  x-axis  for  a  distance  L.  The  top  wire 
starts  at  a  distance  d  on  the  positive  y-axis  and  extends  out  in  the 
positive  x  direction  a  distance  L  to  the  point  x=L,  y=d,  z=0. 

Only  certain  structures  are  capable  of  supporting  a  TEM  mode  of 
propagation.  These  structures  always  have  at  least  two  conductors 
(coaxial  line,  parallel  plate  waveguide,  two  wire  transmission  line). 
One  of  the  most  basic  and  important  assumptions  of  transmission  line 
theory  is  that  only  the  TEM  mode  propagates  down  the  line.  The  re¬ 
strictions  on  frequency  and  wire  separation  are  directly  related  to 
this  assumption.  As  mentioned  in  Chapter  1,  the  electrical  separa¬ 
tion  (physical  separation  as  a  portion  of  a  wavelength)  must  be  small 
enough  so  that  only  the  TEM  mode  propagates  and  all  higher  order  modes 
are  cut  off.  The  TEM  assumption  assures  that  in  any  plane  transverse 
to  the  longitudinal  path  of  the  transmission  lines,  static  electro¬ 
magnetic  field  conditions  are  satisfied.  This  is  very  important  in 
that  it  allows  the  calculation  of  a  per-unit- length  capacitance  and 
a  per-unit -length  inductance.  A  unique  voltage  can  then  be  defined 
between  the  conductors  at  a  point,  x,  anywhere  along  the  line.  Also 
a  unique  current,  equal  but  oppositely  directed  in  each  of  the  two 
conductors  (the  current  in  the  top  wire,  flowing  in.  one  direction, 
is  equal  to  the  current  in  the  bottom  wire,  flowing  in  the  opposite 
direction),  can  be  defined  at  any  point  x,  along  the  line. 
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The  property  of  currents  being  equal  and  opposite  is  referred 
to  as  a  "differential  mode".  (See  FIGURE  2-2 (a) . )  Equal  currents 
that  flow  in  the  same  direction  are  referred  to  as  "common  mode" 
currents  (See  FIGURE  2-2 (b)).  Given  any  arbitrary  pair  of  currents 
flowing  at  the  same  position  x  on  the  line,  these  currents  can  be  de¬ 
composed  algebraically  into  common  and  differential  mode  currents 
(see  FIGURE  2-2 (c)).  The  top  current,  I^x),  is  equal  to  the  sum  of 
the  common  and  differential  mode  currents.  The  bottom  current,  ^(x), 
is  equal  to  the  common  mode  current  minus  the  differential  mode  cur¬ 
rent.  Transmission  line  theory  assumes  and  predicts  only  differential 
mode  currents  due  to  the  TEM  assumption. 

Based  on  the  previous  description,  the  following  is  a  presentation 
of  the  derivation  of  the  transmission  line  equations  and  the  solution 
for  the  terminal  currents.  This  derivation  follows  a  standard  and 
accepted  format  such  as  that  presented  in  [1],  [2],  and  [3].  In  the 
following  development,  the  transmission  line  wires  are  assumed  to  be 
non-insulated ,  perfect  conductors  immersed  in  a  homogeneous,  lossless 
medium. 

Faraday's  law,  in  the  integral,  sinusoidal  steady  state  (time- 
harmonic)  form  provides  the  first  step  in  the  derivation  of  the  trans¬ 
mission  line  equations.  Consider  a  Ax  section  (see  FIGURE  2-3(a)) 
of  the  line.  A  differential  surface  in  the  x-y  plane,  with  the  unit 
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1  I0<*) 

Differential  Mode  Currents 

(a) 


ic<*>» 


Common  Mode  Currents 

(b) 


3 


3 


I,  <x)  =  Ir(x)  +  In(x) 

3 


3 


r2(x)  =  ic(x)— iD(  x) 

(c) 


FIGURE  2-2 


vector,  ds,  normal  to  the  surface,  is  defined  by  the  closed  contour 
t.  The  direction  of  the  contour  is  represented  by  the  unit  vector, 
d£,  which  is  related  to  ds  by  the  right-hand  rule.  The  closed  con¬ 
tour  integration  on  the  left-hand  side  of  Faraday's  law. 


The  surface  integration  on  the  right-hand  side  becomes 
r  -  “  rx+Axr  d 

-jwy  I  H  •  ds  =  -  jaiu  I  H  (x,y)dxdy  (2-3) 

J s  -  x  Jy=0  2 

By  the  assumption  of  perfect  electric  conductors,  the  electric 
field  in  the  x  direction  must  be  zero  at  y=0  and  y=d.  Therefore,  both 
line  integrals  involving  Ex  will  be  zero. 


Voltage  on  the  transmission  line  is  defined  as  (see  FIGURE  2-4) 


=  -f  E  •  dl  =  -f  E  (x,y)dy  =f  E  (x,y)dy  (2-4) 

Jo  Joy  Ja  y 
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FIGURE  2-4 


for  some  fixed  value  of  x.  The  voltage  at  x  +  Ax  becomes 


-  V (x+Ax)  =f  Ey(x+Ax,y)dy  (2-5) 

J  0 


Equation  (2-1)  now  has  the  form 


V  (x+Ax)  -  V(x) 


r  x+Ax  /  d 

=  j<4j)  I  H  (x, 
«'y=0 


y)dxdy 


(2-6) 


Now,  considering  the  left-hand  side  again,  and  dividing  by  Ax 


LIM 

Ax->-0 


N (x+Ax)  -  V(x) 


Ax 


-  k  vw 


(2-7) 


which  gives  the  relationship 


j  ,d 

^  V (x)  =  juy J  Hz(x,y)dy 

»  n 


The  H  field,  or  Hz>  can  be  seen  as  being 
and  a  scattered  field. 


H  =  H1  ♦  HS 
z  z  z 


(2-8) 

the  sum  of  an  incident 

(2-9) 


Now  the  current  on  the  top  wire  can  be  defined  in  terms  of  the  trans¬ 
verse  scattered  H  field.  (See  FIGURE  2-5.) 
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I  (x)  =  f  •  dl 


(2-10) 


Also,  the  per-unit -length  flux  linkage  can  be  defined  in  terms  of 
-*s 


’l'(x)  =  j  •  (ax  X  di2) 


(2-11) 


where  represents  the  direction  of  I(x).  By  definition  [9],  the 
per-unit -length  inductance,  l,  is 

ft  •  (ax  X  “V 
/.*■ 


1  = 


dl. 


Applying  (2-9)  to  equation  (2-8)  we  obtain 


(2-12) 


E  V(x)  = 


/  d  /  d 

j  toy/  H*(x,y)dy  +  jc oul  H^(x,y)dy  (2-13) 


Equation  (2-11)  can  be  written  as 


f d  s 

t|>(x)  *  -  II W  =  v  H  dy 

Jo 


(2-14) 


where  H z  equals  -  Hj,  at  z=0.  Substituting  equation  (2-14)  into  equa¬ 
tion  (2-13)  yields  the  first  transmission  line  equation. 


f  d 

■J  Hx(x, 


fa  V(x)  +  j w£l (x)  =  juv 


y)dy 


(2-15) 
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The  concept  of  continuity  of  current  (conservation  of  charge) 


implies  the  following  relation. 


(2-16) 


Examining  a  differential  volume  (see  FIGURE  2-6(a)),  the  left-hand 
side  of  (2-16)  can  be  divided  into  two  surface  integrations,  with 
normal  surface  vectors  shown  in  FIGURE  2-6(b). 


(2-17) 


The  integral  involving  the  surface  sg  (endcap  surfaces)  represents 
the  axial  conduction  current.  This  current  flows  in  one  end  of  the 
segment  and  out  of  the  other  end;  or  more  properly,  it  represents 
the  currents  flowing  onto  and  off  of  the  segment,  just  off  the  sur¬ 
face  of  the  conductor  (since  perfect  conductors  are  assumed) ,  in  the 
direction  of  the  unit  vector  dsg  (see  FIGURE  2-7 (a)).  This  term  can 
be  represented  as  the  difference  between  the  current  leaving  the  seg¬ 
ment  at  x+Ax  and  the  current  entering  the  segment  at  x. 


I(x+Ax)  -  I(x) 


■/. 


J 


ds 


(2-18) 


The  second  integration  term,  involving  s£  (circumferential  or 
cylindrical  surface)  represents  the  transverse  conduction  current 
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I  (x  +  Ax  ) 


I  (x) 


x  x  +  A 


Conduction  Current 

(a) 


Transverse  Conduction  Current 
(b) 


FIGURE  2-7 
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(see  FIGURE  2-7 (b)). 


(2-19) 


where  o  is  the  conductivity  of  the  medium  and  is  the  E  field  trans¬ 
verse  to  the  conductors.  The  assumption  of  a  lossless  medium  allows 
us  to  neglect  this  term. 


The  integral  on  the  right-hand  side  of  (2-16)  represents  the 
total  charge  enclosed  by  the  differential  volume. 


L 


pdv  =  Q 


(2-20) 


where  p  is  the  charge  density  and  Q  is  the  total  charge  enclosed.  A 
cross  section  of  the  differential  segment  of  transmission  line  is 
shown  in  FIGURE  2-8.  Here,  Q  can  also  be  defined  in  terms  of  the 
scattered  E  field. 


r  x+Ax  * 

=  V  J 

J  x  J  c 


c?T 


,dx 


(2-21) 


where  the  contour,  c,  is  just  off  the  surface  of  the  conductor.  As 
defined  by  equation  (2-4) ,  the  voltage  due  to  the  scattered  field  is 


VS  (x)  =  - 


di. 


(2-22) 
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FIGURE  2-8 
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By  definition  [9],  the  per-unit-length  capacitance  is 


(2-23) 


(2-24) 


With  the  substitution  of  equation  (2-24)  for  the  right-hand  side 
of  (2-16),  we  obtain 

I  (x+Ax)  -  I(x)  -  -  jAx  uc  VS (x)  (2-25) 


The  total  voltage  has  not  yet  been  represented  in  the  above 

5 

equations.  V  (x)  represents  only  the  voltage  due  to  the  scattered 
field.  The  total  voltage  is  the  sum  of  the  voltage  due  to  the  inci¬ 
dent  field  and  the  voltage  due  to  the  scattered  field. 


(2-27) 


VS (x)  =  V(x)  -  V1 (x) 

By  definition  (equation  (2-4)), 

/•d 

-  V1 (x)  =  E*(x,y)dy  (2-28) 

Jo  y 

Substitution  of  (2-28)  and  (2-27)  yields  the  following  from  (2-26). 

^  I(x)  =  -  jwc  V (x)  -  jwc j  Ey(x,y)dy  (2-29) 

This  is  the  second  transmission  line  equation. 

The  two  transmission  line  equations  can  easily  be  put  in  a  form 
suitable  for  solution  by  state  variable  methods. 

J;  V  (x)  =  -  jwil  I(x1  +  Vs(x)  (2-30) 

^  I(x)  =  -  jwc  V(x)  +  Is(x)  (2-31) 

where 


Vs(x)  =  juuj"  H*(x,y)dy 

(2-32) 

r  d  . 

Is(x)  =  -  jwcj  Ey(x,y)dy 

'  A 

(2-33) 

0 
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If  we  define  Z=jw£  and  Y=jtuc  the  following  equation  is  in  state 
variable  form. 


■ 

—  - 

r  *"! 

•  « 

V(x) 

o  -z 

V(X) 

Vs  Cx) 

= 

+ 

joo 

-Y  0 

I(X) 

I  (x) 

Ls  J 

(2-34) 


This  has  the  same  form  as  the  more  general  equation, 


X(t)  =  A  X(t)  +  F(t) 


(2-35) 


often  encountered  in  linear  system  theory.  Here,  X(t)  is  the  state 
vector  whose  entries  are  the  system  states  and  F(t)  is  the  forcing 
vector  whose  entries  are  the  system  inputs. 

The  state  variable  solution  as  presented  in  [4]  is 
A(t-t  )  rt  £(t-x) 

X(t)  =  e~  0  X(t  )  +1  e  F(x)dT  (2-36) 

0  J  to  - 

In  transmission  line  theory, 

A(t)  A 

e  =  4>  Cx)  (2-37) 

where  <p(x)  is  the  chain  parameter  matrix.  It  serves  to  relate  vol¬ 
tages  and  currents  at  different  positions  on  the  line.  For  a  loss- 
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less  transmission  line  [8] 


♦ll(x)  <!>12(x) 

4>2i  (x)  <t>22('X-) 


where 


cos(6x)  -  jvJ l  sin(Bx) 

-  jvc  sin(Bx)  cos(Bx) 


(2-38) 


8  « 


2tt 

X 


and  is  the  propagation  constant  (X  =  wavelength),  v  is  the  velocity 
of  light,  and  l  and  c  are  the  per-unit-length  capacitance  and  induc¬ 
tance  respectively,  as  defined  earlier. 


Specializing  (2-36)  to  the  two  wire  transmission  line  model,  the 
following  relation  is  obtained 

r, 

Yv.(t)l 

dt  (2-39) 


V  (x) 

Yv(0) 

VS(T) 

= 

♦  (X) 

+ 

4>  (x-t) 

I(x) 

L '  II(0)J 

J 

0 

Writing  separate  equations  for  V(x)  and  I(x)  gives: 


V(x)  =  cos(Sx)  V(0)  -  jvH  sin(Sx)  1(0)  +  V'(x) 


(2-40) 


I(x)  =  -  jvc  sin(Bx)  V(0)  +  cos(Bx)  1(0)  +  I'(x)  (2-41) 
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where 


V"(x)  =f  [cos(8(x-x))V  (x)  -  jvZ  sin(6(x-x))I  (x)]dx  (2-42) 
J0 

I'(x)  =J  l-  jvc  sin(6(x-x))Vs(x)  +  cos (6 (x-x)) Ig (x) ]dx  (2-43) 

In  order  to  solve  the  system  of  equations,  (2-40)  and  (2-41), 
the  terminal  conditions  must  be  specified  to  relate  V(0)  and  1(0), 
and  V(x)  and  I(x).  If  we  let  x=L,  the  line  length,  then  V(L)  is  re¬ 
lated  to  I (L)  by  the  load  at  x=L,  Z^. 

V(L)  *  ZLI(l)  (2-44) 

(See  FIGURE  2-9.)  V(0)  is  related  to  1(0)  by  the  load  at  x=0,  ZQ. 

V (0)  =  -  ZQI (0)  (2-45) 

Substituting  (2-44)  and  (2-45)  into  (2-40)  and  (2-41), 


I(L)  =  ~  [-  cos (8L)ZqI (0)  -  jvZ  sin(8L)I(0)  +  V'(L))  (2-46) 


I (L)  =  jvc  sin(8L)ZQI(0)  +  cos(8L)I(0)  +  I'(L)  (2-47) 
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Setting  (2-46)  and  (2-47)  equal  and  solving  for  1(0)  gives  the  fol¬ 
lowing,  which  is  equivalent  to  the  form  presented  in  [1]  specialized 
to  the  two-conductor  case: 


R  "  20  -  ZL  <2-49> 

Also,  the  characteristic  resistance  of  the  line,  R£,  can  be  defined 
as  [8] 

R  =  v  £.  =  — -  (2-50) 

c  vc 


Equation  (2-48)  can  now  be  written  in  the  form 


1(0)  = 


Vg(L)  -  R  I'(L) 


(2-51) 


2 (R  cos(8L))  +  j(R  +g-)  sin(BL) 
C  Rc 


With  the  substitution  of  the  result  of  (2-51)  for  1(0)  (and  R  1(0) 


for  V(0))  into  equation  (2-41),  the  transmission  line  solution  for  the 
differential  mode  current  at  any  point  along  the  line  can  be  obtained, 
including  x=L,  as  given  by  equation  (2-46)  or  (2-47). 
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The  forcing  functions  V'(x)  and  I'(x)  contain  the  terms  V  (x) 
and  I  (x),  given  by  equations  (2-32)  and  (2-33),  which  represent 
sources  induced  by  an  incident  field.  One  common  way  of  supplying 
this  incident  field  is  to  specify  a  uniform  plane  wave  excitation. 

A  uniform  plane  wave  (see  FIGURE  2-10)  has  the  property  that  both 
the  electric  field  vector,  E,  and  the  magnetic  field  vector,  H,  are 
transverse  to  the  direction  of  propagation,  represented  by  the  vector 
?.  These  vectors  are  related  by  the  following  vector  expression. 

E  x  H  =  P  (2-52) 

This  is  TEM  propagation,  as  discussed  earlier. 

The  details  of  how  to  get  from  an  incident  uniform  plane  wave  to 
the  E  and  H  field  components  needed  for  Vs(x)  and  Is(x)  are  provided 
in  Appendix  I.  A  computer  program  was  written  to  generate  transmission 
line  terminal  current  predictions  for  the  two  wire  structure,  when 
excited  by  an  incident  plane  wave.  A  description  of  this  program  is 
given  in  Appendix  II  and  a  listing  is  given  in  Appendix  III. 

Although  the  equations  involved  in  going  from  an  incident  plane 
wave  to  a  prediction  of  the  current  entering  a  load  appear  somewhat 
formidable  and  quite  cumbersome  (and  indeed  they  can  be),  greatly  sim¬ 
plified  forms  are  presented  in  Chapter  4  for  the  special  cases  con- 
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sidered  there.  These  forms  allow  more  insight  into  the  relationship 
between  incident  field  excitation  and  load  currents. 


CHAPTER  3 


THE  METHOD  OF  MOMENTS  FORMULATION 


The  transmission  line  model  provides  a  very  efficient  technique 
for  predicting  the  terminal  currents  induced  by  an  incident  field. 
Although  the  theory  behind  this  model  is  well  known,  the  limits  of 
its  validity  have  not  yet  been  established.  Experimental  data  for 
a  transmission  line  excited  by  an  incident  field  are  very  difficult 
to  obtain.  To  properly  measure  the  induced  currents,  at  relevant 
frequencies,  would  probably  require  the  use  of  an  anechoic  chamber. 
Unfortunately,  this  type  of  facility  is  rarely  available  for  these 
kinds  of  measurements.  However,  in  lieu  of  experimental  data,  a  nu¬ 
merical  solution  technique  called  the  "method  of  moments"  has  been 
used  in  Chapters  4  and  5  to  generate  a  set  of  baseline  data  with  which 
the  transmission  line  model  predictions  can  be  compared. 

The  method  of  moments  has  gained  wide  acceptance  and  credibility 
in  recent  years.  Several  user-oriented  computer  programs  have  been 
developed  to  solve  problems  in  electromagnetics  [10]  (particularly 
in  the  area  of  electromagnetic  compatability) ,  since  Harrington  [11] 
specialized  the  technique  to  electromagnetic  problems.  While  the 
method  of  moments  can  give  very  accurate  results,  the  cost  (computer 
time  and  memory)  of  this  accuracy  can  be  substantial,  as  mentioned  in 
Chapter  1.  This  is  why  computationally  efficient  models,  like  the 
transmission  line  model,  are  very  important. 
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One  reason  that  the  method  of  moments  technique  is  so  expensive 
is  that  very  few  assumptions  are  made.  One  significant  assumption, 
however,  in  both  of  the  formulations  presented  later  in  this  chapter 
is  that  of  "thin  wires".  In  addition  to  the  implications  discussed  in 
earlier  chapters  for  the  thin  wire  assumption,  it  has  an  important 
effect  on  the  enforcement  of  E  field  boundary  conditions.  This  will 
be  discussed  further  in  a  moment. 


In  the  method  of  moments  technique  the  current  distribution  over 
the  entire  structure  is  solved  for  by  enforcing  the  boundary  condition 
that  the  tangential  E  field  should  be  zero  at  any  point  on  the  struc¬ 
ture.  Due  to  the  thin  wire  approximation,  this  condition  is  applied 
only  to  the  axial  component  of  the  E  field  at  the  surface  of  the  con¬ 
ductors.  The  method  of  moments  can  be  thought  of  in  the  following  way. 
The  structure  is  divided  up  into  many  subsections.  The  currents  on 
each  of  the  subsections  can  be  put  into  vector  form. 


I 


1 


I 


(3-1) 
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where  1^  (i=l ,2, . . . ,n)  is  the  current  on  the  i**1  segment  (subsection). 
The  voltage  of  each  segment  can  be  treated  in  the  same  fashion  [10]. 


These  voltages  and  currents  can  then  be  related  by  a  generalized  im¬ 
pedance  matrix,  which  essentially  characterizes  the  structure. 


[V]  =  (Z 1 [ I ] 


(3-3) 


The  values  of  the  entries  in  the  voltage  vector  are  supplied  by  the 
excitation  (incident  field)  for  scattering  problems.  Standard  matrix 
analysis  techniques  can  now  be  used  to  solve  for  the  current  vector. 

[I]  =  [Z]"1  m  (3-4) 

Loading  can  be  accounted  for  by  adding  the  value  of  the  load  to  the 
value  of  the  impedance  for  the  corresponding  segment  on  the  diagonal 
of  the  impedance  matrix.  Losses  in  the  conductors  can  be  handled  in 
the  same  way,  adding  a  load  to  every  segment.  A  description  of  the 


-39- 


method  of  moments  in  terms  of  linear  operators  is  given  in  Chapter  1 
of  [11]. 

Basically,  this  technique  allows  the  complex  integro-differential 
equations  that  are  normally  used  to  characterize  the  current  on  this 
type  of  structure  to  be  put  into  matrix  form  and  solved  as  a  linear 
system.  The  current  on  each  segment  must  be  represented  by  some 
function,  referred  to  as  an  expansion  function.  One  of  the  major  dif¬ 
ferences  in  various  formulations  of  the  method  of  moments  is  the  choice 
of  expansion  functions.  A  second  major  difference  is  the  choice  of 
testing  functions.  The  testing  functions  are  used  to  enforce  the  tan¬ 
gential  E  field  boundary  condition.  The  two  method  of  moments  pro¬ 
grams  to  be  discussed  here  employ  different  expansion  and  testing 
functions.  The  second  program  discussed  uses  a  special  procedure  re¬ 
ferred  to  as  Galerkin's  method  where  the  expansion  and  testing  func¬ 
tions  are  the  same. 

The  first  computer  program  to  be  discussed  is  entitled  WRSMOM, 
and  was  written  by  D.  E.  Warren.  WRSMOM  is  a  program  that  was  de¬ 
signed  specifically  for  configurations  of  straight  thin  wires  [10]. 

It  is  based  on  the  formulation  presented  in  Chapter  4  of  [11].  The 
program  employs  the  "pulse  expansion  and  point  matching"  technique. 

This  refers  to  the  current  expansion  functions  and  the  E  field  testing 
functions  discussed  earlier.  In  this  technique,  the  current  is  assumed 
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to  be  constant  over  each  segment;  or  in  other  words,  the  current  on 
each  segment  is  expanded  as  a  piecewise  constant  pulse.  (See  FIGURE 
3-1 (a).)  The  tangential  E  field  boundary  condition  is  enforced  only 
at  the  "matchpoint"  (hence  "point  matching")  located  at  the  center 
of  each  segment.  The  testing  function  is  therefore  a  form  of  the 
dirac  (or  delta)  function.  The  program  output  includes  the  current 
array,  which  is  the  same  as  the  current  vector  discussed  earlier. 

Each  entry  in  this  vector  gives  the  magnitude  and  phase  of  the  pulse 
representing  the  complex  current  over  the  corresponding  subsection. 

Due  to  the  way  in  which  charge  and  current  are  represented  on 
the  structure,  it  is  important  that  the  segment  length  (or  zone  size) 
remains  fairly  consistent.  Although  segments  of  varying  lengths  can 
be  used,  only  small  changes  in  length  between  adjacent  subsections 
are  acceptable.  To  maintain  the  thin  wire  assumption,  the  zone  size 
should  be  several  times  greater  than  the  radius  of  the  wire  (at  least 
five  times  the  radius  [12]).  Also,  ten  to  forty  subsections  per  wire 
wavelength  are  recommended  for  accuracy  in  [12]. 

The  pulse  expansion  and  point  matching  technique  is  fairly  well 
known  and  popular  because  it  is,  conceptually,  the  most  straightforward 
formulation.  Also,  this  is  the  technique  first  specialized  to  electro¬ 
magnetic  problems  by  Harrington.  However,  other  formulations  using  more 
sophisticated  expansion  and  testing  functions  can  in  some  cases  give 
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(b) 


FIGURE  3-1 
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more  accurate  results  and  better  convergence.  One  such  computer  code 
is  OSMOM  (this  is  not  the  name  given  it  by  the  author,  but  is  one  that 
is  commonly  used)  written  by  J.  H.  Richmond  at  Ohio  State  University 
[13],  [14].  Richmond's  code  has  gained  a  great  deal  of  respect  and 
a  general  acceptance  from  users  of  the  method  of  moments  [10].  In 
fact,  several  individuals  have  used  the  code  as  a  basis  for  modified 
versions.  The  program  is  noted  for  handling  wire  junctions  and  cor¬ 
ners  very  well.  This  makes  it  well  suited  for  application  to  the 
(transmission  line)  structures  considered  in  the  next  few  chapters. 

The  expansion  functions  used  in  OSMOM  are  piecewise  sinusoidal  (some¬ 
times  referred  to  as  sinusoidal  triangles).  (See  FIGURE  3-l(b).) 

These  functions  are  defined  as  follows: 


^(x) 


sin0(x  -  x^  p 

sing  (xi"-'xi_1)  ’  xi-l  -  x  -  xi 


sin0(x^+^  -  x) 

sin8(x.+1  -  x.)  *  xi  -  x  -  xi+l 


(3-5) 


where  8  is  the  propagation  constant  defined  in  Chapter  2.  Here  I^(x) 
is  the  current  approximation  for  an  x-directed  current  (axial  current 
on  an  x-directed  wire),  where  each  expansion  function  spans  two  seg¬ 
ments  and  adjacent  functions  overlap.  The  i**1  entry  in  the  current 
vector  is  the  complex  amplitude  of  the  corresponding  piecewise  sinu¬ 
soidal  expansion  function.  A  wire  which  is  divided  into  n  subsections 
will  have  n+1  entries  in  the  current  vector,  one  entry  foi  the  begin- 
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ning  and  end  of  each  segment  (of  course  the  end  of  one  segment  is  the 
beginning  of  the  next). 

OSMOM  uses  a  special  technique,  referred  to  as  Galerkin's  method, 
so  the  testing  functions  are  the  same  as  the  expansion  functions,  as 
was  mentioned  earlier.  Lumped  loads  can  be  accounted  for  by  adding 
a  diagonal  matrix  of  loads  to  the  generalized  impedance  matrix.  This 
is  the  same  as  with  WRSMOM  except  that  the  loads  are  placed  at  the  end¬ 
points  of  adjoining  segments  (which  correspond  to  the  expansion  func¬ 
tion  peaks)  instead  of  at  the  centers  of  the  segments.  Excitation 
voltages,  which  are  derived  in  the  program  from  the  plane  wave  inci¬ 
dent  field  (as  with  WRSMOM),  are  also  positioned  at  the  peaks  of  the 
sinusoidal  triangles.  Because  of  the  type  of  expansion  function  and 
the  Galerkin  technique,  consistent  zone  size  is  not  as  critical  as 
with  the  pulse  expansion  and  point  matching  procedure.  However,  the 
restrictions  placed  on  subsection  length  by  the  thin  wire  assumption 
still  apply.  For  the  sake  of  comparison,  the  same  zone  structure 
will  be  used  for  both  WRSMOM  and  OSMOM,  applying  the  more  stringent 
requirements  of  the  pulse  expansion  and  point  matching  technique  to 
both,  whenever  possible. 

In  general,  Galerkin's  method  seems  to  give  accurate  results 
with  fewer  subsections  and  better  convergence  than  pulse  expansion 
and  point  matching.  In  addition,  sinusoidal  triangles  seem  to  be 
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very  effective  as  expansion  functions. 


The  transmission  line  structure  presents  some  modeling  problems 
which  particularly  tax  the  method  of  moments  and  especially  the  pulse 
expansion,  point  matching  technique.  With  the  method  of  moments, 
the  terminal  loads  now  become  part  of  the  physical  structure  and  must 
be  modeled  as  wires  containing  lumped  loading.  Thus,  the  structure 
is  made  up  of  two  long  wires  and  two  very  short  wires.  In  order  to 
preserve  any  consistency  in  zone  length  and  at  the  same  time  maintain 
a  reasonable  matrix  size  (a  reasonable  number  of  zones),  very  few 
zones  can  be  used  to  model  the  terminal  wires.  As  suggested  by  [IS], 
in  certain  cases  the  matchpoints  of  the  collocation  [point  matching) 
technique  cannot  be  spaced  sufficiently  dense  to  accurately  sample  the 
E  field.  This  problem  is  compounded  when  the  incident  field  is  such 
that  the  only  excitation  voltages  (in  the  voltage  vector)  are  located 
in  the  loads.  In  this  case  the  highly  localized  driving  field  is  very 
poorly  represented.  In  addition,  if  the  current  magnitude  changes 
rapidly  near  the  termination,  the  pulse  expansion  formulation  cannot 
adequately  model  the  rapidly  changing  subsection  currents.  These 
factors  are  extremely  important  since  we  are  interested  in  an  accurate 
prediction  of  the  currents  entering  the  loads. 

To  summarize,  we  should  expect  the  Galerkin's  method  formulation 
(with  sinusoidal  triangles)  to  give  much  better  accuracy  and  conver- 
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gence  than  the  pulse  expansion,  point  matching  technique  for  trans¬ 
mission  line  structures  with  close  wire  separation,  when  the  same  zone 
configuration  is  used  for  both  solutions.  The  sinusoidal  expansion 
functions  allow  the  tangential  E  field  to  be  more  adequately  repre¬ 
sented  in  and  near  the  terminal  wires  than  do  the  highly  localized 
matchpoints  of  the  collocation  technique.  This  is  especially  true 
in  the  comers  where  the  transmission  lines  meet  the  terminal  wires 
and  it  is  impossible  to  place  a  matchpoint.  For  these  reasons,  the 
transmission  line  results  will  only  be  compared  with  predictions 
made  by  OSMOM  in  Chapter  4.  However,  WRSMOM  will  be  used,  along 
with  OSMOM  in  Chapter  5,  where  several  special  cases  are  considered. 
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CHAPTER  4 


PREDICTION  ACCURACY  OF  THE  TRANSMISSION  LINE  MODEL 


In  the  last  few  chapters  we  have  been  discussing  solution  tech¬ 
niques  for  the  terminal  currents  of  a  general  two  wire  transmission 
line.  The  structure  was  described  by  its  physical  dimensions  of 
length,  wire  separation  and  wire  radius  and  by  the  terminal  loads. 

In  this  chapter,  specific  values  will  be  chosen  for  these  parameters 
and  predictions  made  by  the  different  methods  of  solution  will  be  com¬ 
pared  in  detail.  It  may  be  helpful  to  keep  in  mind  that  the  method 
of  moments  solution  should  establish  a  set  of  baseline  data.  These 
data  have  been  generated  by  a  rigorous  solution  technique  and  are 
accepted  among  experts  in  the  field  as  the  "correct"  solution.  The 
transmission  line  model  provides  a  much  more  efficient  means  of  pre¬ 
diction.  The  effort  here  is  intended  to  validate  the  transmission 
line  results  by  comparison  with  the  baseline  data  set  produced  by  the 
method  of  moments. 

The  dimensions  of  the  structure  used  for  analysis  in  this  chapter 
were  carefully  selected.  It  would  be  an  impossible  task  to  evaluate 
all  configurations  of  transmission  line  networks.  With  the  wide  variety 
of  wire  gauges,  lengths  of  lines,  separations  between  wires,  and  wire 
bundling  techniques,  even  defining  an  "average"  structure  is  quite 
subjective.  However,  several  characteristics  common  to  all  transmission 
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lines  help  guide  the  selection  of  dimensions.  As  was  discussed  in 
Chapter  1,  the  separation  between  the  wires  should  be  many  times  smal¬ 
ler  than  the  overall  length  of  the  line.  For  a  line  with  a  length 
of  one  meter,  a  wire  separation  of  one  centimeter  more  than  satis¬ 
fies  this  requirement.  The  wire  radius  is  assumed  to  be  many  times 
smaller  than  the  separation.  Selecting  the  wire  radius  to  be  one  tenth 
of  a  millimeter  more  than  satisfies  this  requirement  and  also  assures 
the  validity  of  the  thin  wire  assumption.  This  radius  corresponds  ap¬ 
proximately  to  32  gauge  wire,  which  i?  in  the  range  of  gauges  commonly 
used  for  lines  carrying  low  level  signals.  Th_  wire  separation  could 
represent  the  distance  between  two  conductors  in  a  cable  bundle  or  a 
conductor  at  a  height  of  five  millimeters  over  a  ground  plane  (using 
image  theory).  While  this  separation  may  be  wide  for  some  applica¬ 
tions,  remember  that  the  transmission  line  results  will  improve  with 
closer  separations.  Given  the  parameters  of  radius  and  separation, 
the  selected  length  of  one  meter  is  convenient  and  allows  for  the  ob¬ 
servation  of  terminal  current  predictions  with  a  line  that  is  both 
electrically  long  and  electrically  short  over  the  Tange  of  frequencies 
being  investigated.  All  of  the  physical  dimensions  have  now  been 
defined.  The  line  length  is  one  meter;  the  wire  separation  is  one 
centimeter;  and,  the  radius  is  one  tenth  of  one  millimeter  giving  a 
radius  that  is  one  hundred  times  smaller  than  the  separation,  and 
a  separation  that  is  one  hundred  times  smaller  than  the  line  length. 
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The  only  parameters  left  undefined  are  the  load  impedances.  For 
the  sake  of  simplicity,  the  loads  at  each  termination  will  be  equal 
and  resistive.  This  still  leaves  quite  a  range  of  loading  conditions. 
The  choices  can  be  narrowed  by  looking  at  three  cases:  loads  much 
smaller  than  the  characteristic  resistance  of  the  line,  loads  much 
larger  than  the  characteristic  resistance  and  loads  equal  to  the 
characteristic  resistance  (the  latter  are  referred  to  as  matched 
loads).  The  characteristic  resistance  of  the  present  line  under  dis¬ 
cussion  is  552,2262  ohms  as  given  by  equation  (2-50);  so,  this  is 
the  value  to  be  used  for  the  "matched"  case.  The  selection  of  the 
other  two  loading  conditions  is  again,  highly  subjective.  The  choices 
of  50  ohms  for  the  "low  impedance"  case  and  10K  ohms  for  the  "high 
impedance"  case  represent  nominal  input  impedances  for  typical  elec¬ 
tronic  devices.  The  important  point  with  respect  to  these  latter  two 
choices  is  that  one  is  many  times  smaller  than,  and  one  is  many  times 
larger  than,  the  characteristic  resistance  of  the  line.  Other  choices 
in  those  ranges  should  give  results  similar  to  those  obtained  for  the 
selected  values. 

At  this  point  the  structure  has  been  completely  defined  for  solu¬ 
tion  by  the  transmission  line  model.  One  step  remains,  however,  in 
defining  the  structure  for  solution  by  the  method  of  moments.  Each 
wire  must  be  divided  into  subsections,  as  described  in  Chapter  3.  If 
the  terminal  wires  are  characterized  by  one  segment  each,  using  equal 
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length  segments  on  the  remainder  of  the  structure  creates  more  than 
200  total  subsections.  This  would  require  solving  a  system  of  202 
simultaneous  equations  and  unknowns.  This  is  a  very  large  problem 
even  with  the  best  matrix  analysis  techniques.  Numerical  noise  and 
round-off  problems  might  make  the  accuracy  of  the  result  questionable, 
also.  From  a  standpoint  of  economics  (computer  time  and  memory)  and 
numerical  stability,  a  structure  with  SO  total  subsections  is  much 
more  reasonable.  Again  characterizing  the  terminal  wire  as  one  seg¬ 
ment  and  dividing  each  transmission  line  wire  into  24  equal  segments 
gives  a  structure  with  a  total  of  SO  subsections. 

In  Chapter  2  it  was  pointed  out  that  the  transmission  line  model 
solves  for  the  differential  mode  current,  because  of  the  TEM  assump¬ 
tion.  However,  the  method  of  moments  solves  for  the  total  current. 

For  a  valid  comparison  between  the  two  solutions,  the  differential 
mode  current  must  be  extracted  from  the  total  current  given  by  the 
method  of  moments.  This  is  a  simple  algebraic  manipulation.  We  are 
interested  in  the  currents  entering  and  leaving  the  load  at  x=0.  The 
output  from  0SM0M  will  give  a  value  for  the  total  current,  1^,  flowing 
in  the  positive  x  direction  on  the  top  wire  at  x=0  and  a  value  for 
the  total  current,  Ig,  flowing  in  the  positive  x  direction  on  the  bot¬ 
tom  wire  at  x=0.  The  total  currents  can  be  decomposed  into  common 
and  differential  mode  currents  by  the  following  matrix  equation. 
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(4-1) 


The  differential  mode  component,  1^,  obtained  from  the  method  of 
moments  solution  by  this  equation  is  then  compared  to  the  terminal 
current  prediction  made  by  the  transmission  line  model.  The  relation¬ 
ship  between  the  common  mode,  differential  mode  and  the  total  current 
is  discussed  in  Chapter  6. 

Now  that  the  structure  has  been  fully  characterized  for  solution 
by  the  transmission  line  model  and  by  the  method  of  moments,  the  only 
parameters  yet  to  be  defined  are  those  describing  the  incident  field. 

As  mentioned  earlier  and  described  in  more  detail  in  Appendix  I,  the 
excitation  is  provided  by  a  uniform  plane  wave.  Each  of  the  method 
of  moments  programs  and  the  transmission  line  program  will  solve  for 
the  currents  induced  by  a  plane  wave  excitation  with  any  angle  of  in¬ 
cidence  and  any  polarization.  Three  cases  have  been  selected  for 
study.  These  cases  are  representative  of  all  of  the  possible  cases  in 
the  following  ways.  First,  each  coordinate  direction  of  plane  wave 
propagation  is  represented  by  one  of  the  cases.  Secondly,  in  the 
transmission  line  model,  sources  are  induced  by  the  incident  E  field 
and  by  the  incident  H  field  for  an  arbitrary  angle  of  incidence.  The 
cases  being  considered  here  include  an  incidence  angle  and  polarization 
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that  induce  sources  only  from  the  E  field,  that  induce  sources  only 
from  the  H  field  and  one  case  where  both  the  E  field  and  H  field  in¬ 
duce  distributed  sources.  For  each  of  the  following  three  cases  the 
orientation  of  the  incident  plane  wave  will  be  illustrated  and  a  sim¬ 
plified  version  of  the  transmission  line  model  solution  for  1(0), 
specialized  to  that  case,  will  be  derived.  Even  though  the  general 
formulation  of  the  terminal  current  solution  was  used  to  generate  the 
actual  data,  these  specialized  forms  provide  more  insight  into  the 
solution. 

The  first  case  to  be  considered  has  the  plane  wave  incident  on 
the  end  of  the  structure.  (See  FIGURE  4-1.)  The  plane  wave  is  prop¬ 
agating  in  the  positive  x  direction  and  first  strikes  the  transmis¬ 
sion  line  at  the  end  where  x=0;  thus,  it  will  be  referred  to  as  the 
ENDFIRE  case.  The  E  field  is  oriented  in  the  positive  y  direction, 
leaving  the  H  field  oriented  in  the  positive  z  direction.  This  is  the 
only  logical  polarization  (with  the  exception  of  a  negative  y-directed 
E  field  and  a  negative  z-directed  H  field) .  Since  E  and  are  the 
only  field  components  which  induce  currents,  this  polarization  gives 
maximum  pickup  of  both  fields  by  the  structure.  Rotating  the  polari¬ 
zation  90  degrees  either  direction  gives  zero  pickup,  and  therefore  no 
induced  currents.  The  following  is  a  derivation  of  the  terminal  cur¬ 
rent  flowing  out  of  the  load,  onto  the  top  wire,  at  x=0,  for  the  END- 
FIRE  case. 
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Equations  (2-32)  and  (2-33)  provide  the  induced  voltage  and  cur¬ 
rent  sources  respectively  in  terms  of  the  incident  field  components, 

H1  and  E1.  Since  the  total  incident  E  field  is  y-directed  and  the 
z  y 

total  incident  H  field  is  z-directed. 


E 

7 


E  e 
m 


-j  kx 


(4-2) 


H  =  H  e'j^  =  J?  e^kx 
z  m  n 


(4-3) 


where  E  is  the  incident  E  field  magnitude  and  H  ,  the  incident  H 
m  m 

field  magnitude,  is  related  to  Em  by  the  intrinsic  impedance  of  free 
space,  n,  as  discussed  in  Appendix  I.  The  exponential  terms  indicate 
that  the  plane  wave  is  propagating  in  the  positive  x  direction.  From 
equations  (2-32)  and  (2-33)  we  obtain 

V  (x)  =  jwu  /  H*(x,y)dy  =  jwy|  ■-  e"J  x  dy 

s  J  0  J0 

E 

=  jwy  —  d  e"^KX  (4-4) 

rd  rd  -■ k 

I  (X)  =  -  jwcl  Ey (x,y)dy  =  -  jucl  E  e_]  xdy 
s  J0  y  J  0 

=  -  jwc  E  d  e  ^kx  (4-5) 

m 


Substituting  equations  (4-4)  and  (4-5)  into  equation  (2-42)  gives 
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Vg(L)  =  J  [cos(k(L-T))Vs(t)  -  jvZsin(k(L-x))Is(x)]dx 
0 

=  f  [cos (k(L-x))  j  —  E  d  e”^T  (4-6) 

J  n  m 

0 

+  jvilsin(k(L-x)  jwc  Emd  e  ^T]dx 

Evaluation  and  simplification  of  the  right-hand  side  of  (4-6), 

the  details  of  which  are  given  in  Appendix  IV,  yields  the  following 

relation. 


V;(L)  =  j  Emd  sin(kL)  (4-7) 

Substituting  equations  (4-4)  and  (4-5)  into  equation  (2-43)  gives 

r(L)  =  j  [-  jvc  sin(k(L-x))Vs(T)  +  cos(k(L-T))Is(x)]dT 

=  f  -  jvc  sin(k(L-x))  j  —  E  d  e 
J  J  n  m 

0 

-  cos(k(L-x))  juic  Emd  e"-''kT]dx  (4-8) 

Evaluation  and  simplification  of  equation  (4-8)  (see  Appendix  IV  for 
the  details)  yields 

E  d 

(L)  =  ■  j  f  sin(kL)  (4-9) 
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Equations  (4-7)  and  (4-9)  can  now  be  substituted  into  equation  (2-51) 
to  obtain 

j  E  d  sin(kL)  +  j  E  d  ^  sin(kL) 
m  J  m  R 

1(0)  =  - y -  (4-io) 

2(Rcos(kL))  +  j  (R  +  ^-)  sin(kL) 

C  K 

C 

Multiplying  numerator  and  denominator  by  R^  and  rearranging  the  numera¬ 
tor. 


j  E  d(R  +  R)  sin(kL) 

1(0)  =  - - - ^ - = - = -  (4-11) 

2 (R  Rc  cos(kL))  +  j (R^  +  R  )  sin(kL) 

The  imaginary  and  real  terms  in  the  denominator  can  never  be 
zero  at  the  same  frequency,  so  1(0)  is  never  singular.  However, 
nulls  will  occur  at  frequencies  which  cause  the  sine  function  in  the 
numerator  to  be  zero;  sin(kL)  equals  zero  when: 

kL  =  mr,  n  =  0,1,2,...  (4-12) 

As  defined  in  an  earlier  chapter,  k  =  ^,  giving 

-  L  =  mr  (4-13) 

v 

Substituting  2tt  times  frequency,  f,  for  the  radian  frequency,  u>,  and 
rearranging , 
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(4-14) 


c  ,  V  ^ 

f  =  n(2L> 

Substituting  for  the  velocity  of  light  and  line  length  of  one  meter 
gives 


f  =  n(150  MHz)  (4-15) 

This  says  that  zero  nulls  will  occur  at  frequencies  which  are  multi¬ 
ples  of  150  megahertz,  or  when  the  line  length  is  a  multiple  of  a  half¬ 
wavelength. 

At  low  frequencies,  when  k  becomes  small,  the  cosine  function  in 
the  denominator  will  approach  one  and  the  imaginary  term  will  become 
insignificant  as  sin(kL)  approaches  zero,  so  the  denominator  will 
become  constant.  In  the  numerator,  sin(kL)  will  be  approximately  kL. 
Since  the  numerator  will  then  vary  approximately  as  frequency  varies 
and  the  denominator  will  be  essentially  constant,  1(0)  will  vary  ap¬ 
proximately  as  a  linear  function  of  frequency. 

At  high  frequencies,  1(0)  will  reach  a  maximum  value  at  the  mid¬ 
point  between  each  null,  or  when 

kL  =  nir  +  g  (4-16) 
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which  corresponds  to 


f  =  n(150  MHz)  +  75  MHz  (4-17) 

At  these  points  sin(kL)  is  equal  to  one  and  cos(kL)  is  equal  to  zero, 
and  equation  (4-11)  becomes 


E  d(R  +  R) 

1(0)  =  ±.c  .  (4-18) 

(  c  +  R  > 

giving  the  maximum  current  value,  which  is  a  constant. 

Anticipating  the  behavior  just  described,  the  transmission  line 
program  listed  in  Appendix  III  was  used  to  generate  terminal  current 
predictions  foT  the  three  loading  conditions  described  earlier.  The 
method  of  moments  solution  was  produced  by  the  program  OSMOM.  Trans¬ 
mission  line  model  predictions  are  represented  by  the  solid  curve  and 
OSMOM  predictions  are  represented  by  the  "0's"  at  discrete  points  on 
all  plots  in  this  chapter. 

PLOTS  4-1 (a)  and  4-1  (b)  show  the  magnitude  and  phase  respectively 
for  terminal  currents  at  x=0  with  matched  loads  and  ENDFIRE  excitation. 
The  transmission  line  model  results  compare  quite  favorably  with  the 
method  of  moments  results  over  the  fairly  wide  frequency  range  shown, 

10  megahertz  to  1  gigahertz.  PLOTS  4-2 (a)  and  4-2 (b)  show  the  same 
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configuration  over  the  lower  frequency  range  of  100  kilohertz  to  10 
megahertz.  Two  important  characteristics  should  be  pointed  out  here. 
First,  the  linear  response  of  the  transmission  line  model  predictions 
occurs  as  was  anticipated.  Secondly,  the  method  of  moments  technique 
breaks  down  numerically,  as  expected.  Predictions  by  the  method  of 
moments  could  not  be  obtained  for  frequencies  below  those  represented 
and  those  shown  are  not  particularly  good. 

PLOTS  4-3 (a)  and  4-3 (b)  show  the  magnitude  and  phase  respectively 
for  the  same  conditions  as  PLOTS  4-1 (a)  and  4-1 (b)  except  that  here  the 
low  impedance  loads  are  used  instead  of  matched  loads.  The  transmis¬ 
sion  line  model  predictions  are  generally  good.  However,  some  slight 
discrepancies  occur  between  the  method  of  moments  results  and  the 
transmission  line  model  results  as  frequency  increases.  These  dif¬ 
ferences  might  be  caused  by  reradiation  off  of  the  terminal  wires,  ac¬ 
counted  for  in  the  method  of  moments  solution.  That  is,  the  method  of 
moments  models  the  loads  as  part  of  the  physical  line,  as  discussed 
in  Chapter  3.  For  low  impedance  loads,  relatively  large  currents 
flowing  on  the  terminal  wires  radiate  as  antennas.  This  reradiation 
by  the  terminal  wires  is  accounted  for  in  the  method  of  moments  solu¬ 
tion,  but  not  in  the  transmission  line  model.  This  effect  is  similar 
to  that  encountered  in  [16].  It  was  noted  in  [16]  that  transmission 
line  theory  did  not  predict  resonances  for  currents  induced  in  the 
terminations  of  a  transmission  line  where  the  loads  were  replaced  by 
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short  circuits  at  both  ends.  Further,  it  was  inferred  that  these 
resonances  do  exist  for  the  short  circuit  case  (which  is  an  extreme 
of  the  low  impedance  case)  as  predicted  by  other  solution  methods 
(such  as  the  method  of  moments).  This  explanation  is  enforced  by  the 
predictions  for  higher  impedance  loads,  where  discrepancy  does  not 
occur.  High  impedance  loads  would  tend  not  to  radiate  as  freely  as 
low  impedance  loads.  As  the  impedance  of  the  loads  is  increased, 
eventually  the  termination  will  begin  to  look  like  an  open  circuit 
and  the  transmission  line  will  behave  as  two  separate  wires  in  free 
space. 

PLOTS  4-4 (a)  and  4-4 (b)  are  the  low  frequency  predictions  for 
the  same  low  impedance  configuration  of  PLOTS  4-3 (a)  and  4-3  (b).  The 
numerical  solution  of  Method  of  Moments  is  inaccurate  in  this  range.  PLOTS 
4-5 (a  and  b)  and  4-6 (a  and  b)  show  the  predictions  for  high  impedance 
loads.  The  results  are  very  good  as  was  with  the  matched  case.  The 
plots  for  the  ENDFIRE  case  are  summarized  in  TABLE  4-1 (a). 

The  second  case  to  be  considered  has  the  plane  wave  incident  from 
the  side  of  the  structure  (see  FIGURE  4-2),  this  it  will  be  referred 
to  as  the  SIDEFIRE  case.  The  wave  is  propagating  in  the  positive  y 
direction  striking  the  bottom  wire  first.  From  the  transmission  line 
model  standpoint,  the  E  field  cannot  contribute  any  induced  sources  be¬ 
cause  it  does  not  have  a  y  component.  However,  choosing  an  x-directed 
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FIGURE  4-2 


E  field  aligns  the  H  field  with  the  negative  z  axis,  giving  it  maxi¬ 
mum  pickup  by  the  structure.  As  with  the  ENDFIRE  case,  rotating  the 
polarization  90  degrees  either  direction  would  cause  no  induced 
sources  to  be  produced  in  the  transmission  line  model  and  hence,  no 
induced  currents. 


For  the  SIDEFIRE  case  E1  and  H1  are  as  follows. 

y  z 


E1 

y 


0 


H 


i 

z 


E 

jn 

n 


e 


-Jky 


(4-19) 


(4-20) 


Substituting  these  values  into  equations  (2-32)  and  (2-33), 


Vs(x)  =  jwu 


E  -v 

j?e'Jky 

n 


dy 


nk 


E 

m 


(e'jkd 


1) 


(4-21) 


and 


Is(x)  =  0 

Equation  (4-21)  can  be  simplified  by  noting  that 


(4-22) 
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•<\v. 


?,/***+  -- 


and 


k  =  u  t/TTe 


so  that 


_  _  MM  _  1 

nk  fi  O  4TF ) 

Equation  (4-21)  now  becomes 

V  =  E  (e”^kd  -  1)  (4-23) 

S  m 

With  the  substitution  of  equations  (4-23)  and  (4-22),  equation  (2-42) 
becomes 


r(L)  =  J  cos(k(L-x))Em(e"jkd  -  l)dx 
0 

=  y  (e"jkd  -  l)sin(kL)  (4-24) 

Again  substituting  equations  (4-23)  and  (4-22),  equation  (2-43)  be¬ 
comes 


I'(L)  ■  f  -  jvc  sin(k(L-t))Em(e"^kd  -  l)dx 

JQ 

-jE  ... 

=  jr-g2  (e"JKG  -  1)(1  -  cos(kL))  (4-25) 
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Equations  (4-24)  and  (4-25)  can  now  be  substituted  into  equation  (2-51) 

to  yield  the  following  relation. 

E  .  E  R  . 

y  (e"ikd-l)sin(kL)  +  j  jj2_  (e'jkd-l)  (1-cos (kL)) 

1(0)  -  - y1 - 

2(R  cos(kL))  +  j(Rc  +  5_)  sin(kL) 


Y  (e‘ikd-l)[Rc  sin(kL)  +  jR(l-cos(kL))] 
2  R  Rc  cos (kL)  +  j (R^  ♦  R2)  sin(kL) 


(4-26) 


Again,  the  denoroin.' .:t  of  equation  (4-26)  can  never  be  zero,  as  was 
the  case  with  equation  (4-11).  Nulls  will  occur  when  sin(kL)  equals 
zero  and  cos(kL)  equals  one,  causing  the  numerator  to  become  zero. 
This  condition  corresponds  to 


kl  =  n(2ir)  ,  n=0,l,2, . . .  (4-27) 


or 


f  *  n(300  MHz)  (4-28) 

which  indicates  that  nulls  occur  at  frequencies  which  are  multiples 
of  300  megahertz,  or  when  the  line  length  is  a  multiple  of  one  wave¬ 
length.  The  value  of  the  current  halfway  between  each  null  is  not 
constant,  as  it  had  been  with  the  ENDFIRE  case.  These  current  values 


occur  when 


kL  *  n(2ir)  +  7t 


(4-29) 


or 


£  =  n(300  MHz)  +  150  MHz 


(4-30) 


At  these  frequencies,  cos(kL)  equals  minus  one  and  sin(kL)  equals 
zero,  yielding  the  following  from  equation  (4-26). 


1(0) 


E„Ce-Jk,i  -  1) 

TF 


(4-31) 


This  equation  yields  the  current  halfway  between  any  two  nulls  and 

e’jkd  -  1 

varies  with  frequency  as  - ^ - . 

At  low  frequencies,  the  imaginary  term  will  become  insignificant  com¬ 
pared  to  the  real  term  in  the  denominator  of  equation  (4-26) .  In  the 
numerator,  cos(kL)  will  approach  one,  sin(kL)  will  become  kL,  and 
e~^kd  will  become  1  -  jsin(kd)  or  1  -  jkd.  Therefore  equation  (4-26) 
will  become 


1(0)  = 


-  j  E  kdL 
m 

25 


(4-32) 
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which  suggests  that  current  will  become  approximately  a  linear  func¬ 
tion  of  frequency  at  low  frequencies,  as  it  did  in  the  ENDFIRE  case. 

The  plots  of  data  for  the  SIDEFIRE  case  follow  the  same  format 
as  the  plots  for  the  ENDFIRE  case.  Three  loading  conditions  were  again 
examined  for  two  different  ranges  of  frequency.  The  plots  are  iden¬ 
tified  in  TABLE  4-1  (b) .  Each  plot  number  has  two  corresponding  plots, 
(a)  and  (b) ,  showing  magnitude  and  phase,  respectively. 

The  matched  case  shows  very  good  correlation  between  the  method 
of  moments  and  the  transmission  line  model  results.  Reradiation  is 
again  observed  in  the  low  impedance  case  causing  some  discrepancies 
at  higher  frequencies.  The  current  peaks  observed  for  the  high  im¬ 
pedance  case  are  not  infinite  in  magnitude.  These  currents  adhere  to 
the  values  predicted  by  equation  (4-31) .  The  low  frequency  results 
again  illustrate  the  numerical  breakdown  of  the  method  of  moments. 

That  is,  numerical  noise  and  round-off  error  in  the  computer  begin 
to  affect  the  solution  significantly. 

The  third  case  to  be  considered  has  the  plane  wave  arriving  at 
all  parts  of  the  transmission  line  simultaneously.  The  wave  is  prop¬ 
agating  in  the  negative  z  direction.  The  phase  front  of  the  wave 
strikes  the  structure  broadside;  hence,  this  will  be  referred  to  as 
the  BROADSIDE  case.  (See  FIGURE  4-3.)  Again  the  polarization  is  such 
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that  the  E  field  is  oriented  for  maximum  pickup  by  the  structure. 
Choosing  a  y-directed  E  field  leaves  an  x-directed  H  field.  However, 
the  H  field  contributes  no  induced  sources  since  it  cannot  have  a 
z  component.  This  configuration  yields  the  following. 

E  =  E  e^kz  =  E  (4-33) 

y  m  m  1  ' 

Hz  =  0  (4-34) 

i  kz 

In  equation  (4-33) ,  eJ  is  always  equal  to  one  because  the  structure 
is  located  at  z=0. 

Substituting  these  values  into  equations  (2-32)  and  (2-33) 
yields 


Vs(x)  =  0 


(4-35) 


-  xs 


rd 

=  -  JU)C  I 
J  0 


E  dy  =  -  jujc  E  d 
m  J  J  m 


(4-36) 


With  the  substitution  of  equations  (4-35)  and  (4-36)  into  (2-42)  we 
obtain  the  following 
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V'(L)  =  -  view  E  d 
s  m 


sin(k(L-r))dx 


=  -  E  d  (I  -  cosCkL)) 
m 


(4-37) 


where  the  relation  view  =  k  was  used,  as  illustrated  in  Appendix  IV. 
Substituting  equation  (4-35)  and  (4-36)  into  equation  (2-43)  yields 


I'(L)  =  -  jioc  Emd  f  cos(k(L-T))dx 
J0 


E  d 

•  -  j  -pr-  sin(kL) 
c 


(4-38) 


where  the  relation  —  =  R  was  used  as  was  illustrated  in  Appendix  IV. 

uc  c  rr 


Equations  (4-37)  and  (4-38)  can  now  be  substituted  into  equation 
(2-51)  to  obtain  the  following  equation  for  1(0). 

-  E  d(l-cos  (kL))  +  j  Ed  sin(kL) 
in  m  k 

1(0)  =  - - 

2  (Rcos(kL))  +  j (R  +  §-)  sin(kL) 

C  K 

c 

E  d(j  R  sin(kL)  -  R  (1  -  cos(kL))] 

=  J2 - ^ - 5 -  (4-39) 

2  R  Rc  cos(kL)  +  j  (R^  +  R  )  sin(kL) 


As  with  both  of  the  previous  cases,  the  denominator  of  equation 
(4-39)  can  never  be  zero  so  the  equation  will  never  be  singular.  1(0) 
will  however  be  zero  whenever 
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kL  *  n2ir 


(4-40) 


0,1,2,... 


or 


f  =  n(300  MHz)  (4-41) 

At  frequencies  that  are  multiples  of  300  megahertz  the  line  length 
is  a  multiple  of  one  wavelength.  Current  values  at  intervals  half 
way  between  nulls  occur  when 


kL  =  n2w  +  it 


(4-42) 


or 


f  =  n(300  MHz)  +  150  MHz  (4-43) 

At  each  of  these  frequencies  the  current  will  have  the  constant  value 
given  by  the  following  equation,  which  was  obtained  by  evaluating 
equation  (4-39)  at  the  frequencies  given  by  equation  (4-43) . 

E  d 

1(0)  =  -J-  (4-44) 

The  linear  behavior  observed  for  the  ENDFIRE  and  SIDEFIRE  cases 
at  low  frequencies  occurs  here  for  the  BROADSIDE  case  also.  At  low 
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frequencies  equation  (4-39)  becomes 


1(0) 


jE  dkL 
J  m 

'  2R  “ 


(4-45) 


This  equation  was  obtained  by  the  same  method  used  to  obtain  equation 
(4-32).  It  shows  that  the  linear  relationship  between  1(0)  and  fre¬ 
quency  occurs  here  in  the  BROADSIDE  case  as  it  did  in  the  two  previous 
cases . 


Plots  of  the  data  for  the  BROADSIDE  case  with  the  three  loading 
conditions  and  two  ranges  of  frequency  are  identified  in  TABLE  4-1 (c). 
The  matched  case  again  shows  excellent  correlation  between  the  method 
of  moments  and  the  transmission  line  model  solutions.  The  effects  of 
reradiation  are  not  as  noticeable  as  they  were  in  the  previous  cases 
for  low  impedance  loading.  This  is  probably  related  to  the  fact  that 
the  excitation  over  the  entire  structure  is  at  the  same  phase.  The 
peaks  observed  for  the  low  impedance  case  are  similar  to  those  observed 
for  the  high  impedance  loading  in  the  SIDEFIRE  case.  However,  here 
each  peak  has  the  same  value,  given  by  equation  (4-44).  The  results 
for  high  impedance  loading  for  the  BROADSIDE  case  are  very  good.  However, 
at  low  frequencies  and  for  all  loading  conditions,  the  Method  of  Moments 
predictions  appear  to  be  inaccurate. 
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In  Chapter  5,  the  structure  used  in  this  chapter  is  modified 
to  illustrate  several  important  points.  All  of  the  plots  in  this 
chapter  compare  only  differential  mode  currents.  In  Chapter  6,  the 
common  mode  currents  and  their  relation  to  the  prediction  of  total 
currents  is  discussed. 
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CHAPTER  5 


SPECIAL  CONSIDERATIONS  FOR  MODIFIED  STRUCTURES 

The  transmission  line  model  predictions  for  the  structure  con¬ 
sidered  in  Chapter  4  showed  a  slight  deviation  from  the  method  of 
moments  solution,  when  the  line  was  terminated  with  low  impedance 
loads.  Reradiation  from  the  terminal  wires  was  discussed  as  a  pos¬ 
sible  cause  for  the  differences  between  the  two  solutions.  If  re¬ 
radiation  is  causing  this  discrepancy,  its  effects  should  become  more 
noticeable  as  the  wire  separation  is  increased.  The  longer  the  termi¬ 
nal  wire,  the  greater  the  coupling  between  this  wire  and  the  incident 
field;  hence,  the  induced  current  is  larger,  causing  greater  reradia¬ 
tion. 


To  illustrate  this  point,  consider  the  SIDEFIRE  case  of  Chapter 
4,  with  low  impedance  loading.  The  solution  of  this  case  for  a  wire 
separation  of  0.01  meter  is  shown  by  PLOT  4-9.  Here  the  effects  of 
reradiation  near  the  nulls  are  noticed.  These  effects  are  greatly 
magnified  when  the  wire  separation  is  increased  to  0.05  meter,  as 
shown  by  PLOT  5-1.  Here  the  predictions  by  the  transmission  line 
model  and  the  method  of  moments  differ  significantly. 

Although  most  of  the  noticeable  differences  between  the  two  solu¬ 
tions  of  PLOT  5-1  are  probably  caused  by  reradiation,  increased  separa- 
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PLOT  5-1  LINE  DIMENSIONS  (METERS) 

-  LENGTH  =  1.0 


tion  creates  other  problems.  As  the  transmission  line  wires  are 
moved  further  apart,  gradually  the  transmission  line  assumptions  (in 
particular  the  TEM  assumption)  begin  to  break  down.  To  examine  these 
problems,  however,  the  effects  of  reradiation  must  be  removed.  To 
accomplish  this,  the  load  is  replaced  by  an  open  circuit.  This  entire¬ 
ly  eliminates  terminal  wire  reradiation  by  eliminating  the  terminal 
wire.  The  structure  is  now  modeled  by  the  method  of  moments  as  two 
parallel  wires.  Obviously,  the  terminal  currents  can  no  longer  be 
compared  because  they  will  be  zero.  Therefore,  the  differential  mode 
currents  at  equally  spaced  points  along  the  transmission  line  will  be 
compared.  Again  using  24  subsections  on  each  wire,  the  method  of 
moments  solution  (OSMOM)  provides  the  value  of  the  current  at  23 
points,  interior  to  the  terminals,  on  each  wire  (25  points,  including 
the  zero  current  value  obtained  at  the  end  of  each  wire).  The  dif¬ 
ferential  mode  current  can  be  extracted  from  the  total  current  solu¬ 
tion,  as  was  done  in  Chapter  4.  The  transmission  line  model  predic¬ 
tions  can  be  obtained  at  these  points  by  the  use  of  the  chain  para¬ 
meter  matrix.  The  derivation  of  the  transmission  line  model  solu¬ 
tion  of  I(x)  for  the  SIDEFIRE  case  with  open  circuit  terminations 
follows.  The  SIDEFIRE  case  was  selected  because  the  ENDFIRE  and 
BROADSIDE  cases  produce  zero  differential  mode  currents  for  all  x. 

This  is  shown  in  Appendix  V. 

From  equation  (2-41),  with  x=L,  we  obtain 
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I(L)  =  -  jvc  sin(gL)V(0)  +  cos(BL)I(0)  +  I'(L)  (5-1) 


Because  the  termination  is  an  open  circuit. 


1(0)  =  I (L)  =  0 


(5-2) 


Therefore,  equation  (5-1)  becomes 


0  =  -  jvc  sin(8L)V(0)  +  I'(L) 


(5-3) 


or 


i;(L) 


(5-4) 


Equation  (2-41),  for  a  general  value  of  x,  gives 


I(x)  =  -  jvc  sin(8x)V(0)  +  I'(x)  (5-5) 


Substituting  equation  (5-4)  into  (5-5)  and  simplifying  yields 


«» -  •  Inifc  ’;<L>  ♦  vsM 


(5-6) 


The  values  for  I'(L)  and  I'(x)  can  be  obtained  from  equation  (4-25), 
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-j  E  ... 

I'(L)  =  TIT  (e"J  -  Dd  -  cos(kL))  (5-7) 

5  c 

-j  E  ... 

rw  =  Yir  (e"J  -  !)(!  -  cos(kx))  (5-8) 
c 

Equation  (5-6)  gives  the  value  of  the  differential  mode  current  at 
any  value  of  x  along  the  line. 

The  transmission  line  model  and  the  method  of  moments  provide 
solutions  for  currents  on  the  open  circuit  structure  at  wire  separa¬ 
tions  of  0.01,  0.05,  and  0.10  meter.  Current  magnitudes  versus  posi¬ 
tion  on  the  line  for  the  three  separations  are  shown  by  PLOTS  5-2, 

5-3,  and  5-4.  TABLE  5-1  lists  the  percent  error  in  the  transmission 
line  model  solution  for  each  point  where  the  method  of  moments  pro¬ 
vides  a  current  value.  The  error  is  very  small  for  the  closest  se¬ 
paration  (this  is  the  same  structure  considered  in  Chapter  4,  ex¬ 
cept  with  an  infinite  impedance  load)  but  gets  progressively  worse 
as  separation  increases.  This  clearly  indicates  that  reradiation  is 
not  the  only  problem  brought  on  by  wide  separations  between  the 
wires.  The  conditions  necessary  for  the  TEM  assumption  to  be  valid 
are  apparently  not  being  satisfied,  for  these  wider  separations. 

Up  to  this  point,  the  method  of  moments  solution  has  been  sup¬ 
plied  solely  by  the  computer  program  0SM0M.  In  Chapter  3  the  restric¬ 
tions  of  the  pulse  expansion  and  point  matching  technique  used  in 
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PLOT  5~2  LINE  DIMENSIONS  (METERS) 

-  LENGTH  -  1.0 

EXCITATION:  SIOEFIRE  SEPARATION  =  0.01 

FREQUENCY®  525  MHZ  RADIUS  =  0.0001 
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PLOT  5-4  LINE  DIMENSIONS  (METERS) 

-  LENGTH  »  1.0 

EXCITATION:  SIOEFIRE  SEPARATION  =>  0.10 

FREQUENCY-  525  HHZ  RADIUS  -  0.0001 
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°0.  00  0.14  0.29  0.43  0.57  0.71  0.86  1.00 

POSITION  ON  WIRE  (X) 


PERCENT  ERROR 


X 

d=0.01 

d=O.OS 

d=0. 1 

0.00000 

0.00 

0.00 

0.00 

0.04167 

-0.05 

1.45 

4.54 

0.08333 

0.82 

3.69 

7.08 

0.12500 

1.06 

4.34 

7.83 

0.16667 

1.15 

4.58 

8.09 

0.20833 

1.18 

4.64 

8.10 

0.2S000 

1.16 

4.56 

7.90 

0.29167 

1.09 

4.28 

7.40 

0.33333 

0.90 

3.61 

6.25 

0.37500 

0.34 

1.68 

3.06 

0.41667 

-4.27 

-12.48 

-18.26 

0.45833 

5.03 

19.49 

35.31 

0.50000 

3.75 

14.25 

25.10 

0.54167 

5.01 

19.48 

35.30 

0.58333 

-4.21 

-12.47 

-18.25 

0.62500 

0.35 

1.68 

3.06 

0.66667 

0.91 

3.61 

6.25 

0.70833 

1.10 

4.28 

7.40 

0.75000 

1.17 

4.55 

7.90 

0.79167 

1.18 

4.64 

8.10 

TABLE  5-1 
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PERCENT  ERROR 

X 

d=0.01 

d=0. 05 

d-0.1 

0.83333 

1.15 

4.58 

8.08 

0.87500 

1.05 

4.33 

7.82 

0.91667 

1.81 

3.68 

7.07 

0.95833 

-0.05 

1.44 

4.53 

1.00000 

0.00 

0.00 

0.00 

TABLE  5-1  (con't.) 


WRSMOM,  consistant  zone  size  in  particular,  were  discussed.  The  struc¬ 
ture  used  in  Chapter  4  involved  a  separation  between  wires  that  would 
have  required  a  total  of  at  least  202  subsections  if  all  zones  were 
of  equal  length.  This  creates  a  very  large  matrix  system  as  men¬ 
tioned  earlier.  However,  with  the  introduction  of  a  structure  invol¬ 
ving  wider  wire  separation,  solution  by  WRSMOM  is  feasible. 

For  reference  sake,  the  three  cases  of  excitation  (ENDFIRE, 
SIDEFIRE,  and  BROADSIDE)  are  solved  for  matched  loads  by  WRSMOM, 
for  the  same  structure  used  in  Chapter  4.  This  structure  had  a  sep¬ 
aration  of  0.01  meter  and  the  method  of  moments  model  used  24  zones 
on  each  line  and  1  zone  for  each  terminal  wire.  It  must  be  pointed 
out  that  the  requirements  necessary  for  a  valid  solution  by  the  point 
matching,  pulse  expansion  technique  (used  in  WRSMOM)  have  not  been 
met.  These  predictions  serve  only  as  reference  for  results  yet  to 
be  discussed.  The  current  magnitudes  predicted  by  WRSMOM  have  been 
added  to  PLOTS  4-1,  4-7,  and  4-13  and  are  shown  in  PLOTS  5-5,  5-6, 
and  5-7,  respectively.  The  solution  by  WRSMOM  is  inaccurate  as  ex¬ 
pected,  for  the  reasons  discussed  at  the  end  of  Chapter  3. 

The  structure  is  now  changed  by  increasing  the  wire  separation 
to  0.05  meter.  Dividing  each  transmission  line  wire  into  20  segments 
and  modeling  each  terminal  wire  as  one  segment  characterizes  the 
structure  with  42  equal  length  zones  for  solution  by  the  method  of 
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PLOT  5  -5  LINE  DIMENSIONS  (HETERS) 

-  LENGTH  »  1.0 

EXCITATION*  ENDFIRE  SEPARATION  »  0.01 

LOADINGS  HATCHEO  RADIUS  -  0.0001 
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FREQUENCY 


FREQUENCY 


moments.  Removing  the  abrupt  change  in  subsection  length  found  in 
the  previous  structure  improves  the  WRSMOM  results  considerably. 

These  results,  for  the  ENDFIRE,  SIDEFIRE,  and  BROADSIDE  cases,  are 
shown  by  PLOTS  5-8,  5-9,  and  5-10  respectively.  The  two  method  of 
moments  solutions  compare  quite  favorably,  even  though  the  E  field  is 
only  being  "matched"  at  one  point  in  the  load  with  the  WRSMOM  solu¬ 
tion.  The  OSMOM  solution  should  be  more  accurate  for  an  additional 
reason.  In  OSMOM,  the  elements  in  the  impedance  matrix  are  calculated 
using  the  rigorous  closed-form  impedance  expressions  in  terms  of 
exponential  integrals.  This  method  of  calculation  provides  a  much 
more  accurate  (but  expensive)  solution  than  the  numerical  integra¬ 
tion  techniques  (such  as  Simpson's  integration  rule)  which  are  nor¬ 
mally  used  to  generate  the  entries  in  the  matrix.  Both  method  of 
moments  solutions  are  sensitive  at  frequencies  near  the  current  nulls, 
because  of  the  rapidly  changing  currents.  It  is  not  surprising 
that  predictions  differ  here.  In  general,  the  predictions  from  all 
three  solution  techniques  compare  well  for  all  three  cases  of  exci¬ 
tation. 
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PLOT  5-9  LINE  DIMENSIONS  (HETERS) 

LENGTH  =  1.0 

EXCITATION :  SIOEFIRE  SEPARATION  =  0.05 

LOADING!  HATCHED  RADIUS  ■  0.0001 
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CHAPTER  6 


DIFFERENTIAL  MODE  AND  COMMON  MODE  CURRENT  COMPARISONS 


The  terms  "common  mode"  and  "differential  mode"  were  defined  in 
Chapter  2,  and  have  been  used  repeatedly  since  then  when  discussing 
currents  flowing  on  the  structure.  The  sum  or  difference  of  these 
currents  gives  the  total  current  flowing  on  the  top  or  bottom  wire 
at  a  specified  value  of  x,  as  defined  by  the  two  equations: 


I 


T 


rC  +  ID 


(6-1) 


I 


B 


(6-2) 


These  currents  are  illustrated  in  FIGURE  2-2.  The  common  and  differen 
tial  mode  currents,  in  terms  of  the  top  and  bottom  currents,  are 
given  by  equation  (4-1). 

The  transmission  line  model,  by  its  assumption  of  TEM  propaga¬ 
tion,  solves  for  the  differential  mode  currents  only.  This  was  dis¬ 
cussed  in  detail,  in  Chapter  2.  The  method  of  moments  formulation 
assumes  only  that  the  currents  are  axially  directed,  so  it  gives  a 
total  current  solution.  This  solution  is  decomposed  into  common  and 
differential  mode  currents  for  comparison  to  the  transmission  line 
model  solution,  as  described  in  Chapter  4.  In  this  chapter,  the  rel- 
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evance  of  the  differential  mode  current  prediction  is  discussed. 

First  it  might  be  helpful  to  gain  some  insight  into  the  common 
mode.  The  common  mode  current,  sometimes  referred  to  as  the  antenna 
mode  current,  can  be  though  of  as  the  "net"  current  passing  through  a 
plane  transverse  to  the  transmission  line.  This  is  why  there  are  no 
common  mode  currents  in  the  transmission  line  model;  the  sum  of  the 
currents  passing  through  any  transverse  plan  is  zero  (or,  zero  net 
current).  We  should  expect  the  common  mode  currents  to  be  larger 
for  the  SIDEFIRE  case,  than  for  the  BROADSIDE  or  ENDFIRE  case,  due 
to  the  orientation  of  the  E  field.  In  the  SIDEFIRE  case  the  E  field 
is  aligned  (parallel)  with  the  transmission  line  wires,  allowing  for 
maximum  pickup.  In  the  other  two  cases  the  E  field  is  perpendicular 
to  the  transmission  line  wires.  (It  should  be  pointed  out  that  this 
orientation  is  exactly  opposite  from  the  E  field  orientation  needed 
to  induce  sources  in  the  transmission  line  model.  There,  maximum  E 
field  pickup  was  achieved  with  the  E  field  orientated  transverse  to 
the  wires,  as  in  the  ENDFIRE  and  BROADSIDE  cases.)  The  results  from 
all  of  the  cases  discussed  in  Chapters  4  and  5  confirm  this  expecta¬ 
tion.  Common  mode  currents  for  the  ENDFIRE  and  BROADSIDE  cases  were 
several  orders  of  magnitude  smaller  than  the  differential  mode  cur¬ 
rents  in  solutions  by  the  method  of  moments  codes.  They,  for  the 
most  part,  are  probably  a  result  of  numerical  "noise"  generated 
during  the  matrix  computations.  However,  for  the  SIDEFIRE  case,  the 
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common  mode  currents  were  the  same  order  of  magnitude  as  the  differen¬ 
tial  mode  currents.  In  [16]  it  is  pointed  out  that  the  relationship 
between  common  and  differential  modes  currents  depends  upon  the  ter¬ 
minations  and  may  be  a  function  of  frequency.  From  the  results  just 
discussed,  it  would  seem  that  the  incident  field  excitation  is  also 
very  important. 

This  raises  a  very  important  issue.  If  the  common  mode  currents 
are  as  large  as  the  differential  mode  currents,  of  what  value  are  the 
transmission  line  model  predictions?  For  an  arbitrary  value  of  x 
along  the  line,  the  transmission  line  model  solves  only  for  the  dif¬ 
ferential  mode  current,  as  stated  several  times.  The  total  current 
on  the  line  is  therefore  not  necessarily  represented.  However,  trans¬ 
mission  line  theory  views  the  terminal  loads  as  lumped  parameters. 

So,  at  the  terminations,  the  common  mode  currents  on  the  top  and  bot¬ 
tom  wires  cancel.  Therefore,  the  net  terminal  current  passing  through 
each  load  is  the  differential  mode  current. 

This  says  that  the  transmission  line  model  should  give  an  accurate 
solution  for  the  total  load  current.  This  last  statement  should  be 
qualified  with  two  restrictions  which  have  been  pointed  out  in  pre¬ 
vious  chapters.  At  high  frequencies,  it  is  no  longer  valid  to  assume 
that  the  load  will  behave  as  a  lumped  parameter.  Wide  separations  be¬ 
tween  the  wires  will  also  cause  this  assumption  to  be  questionable. 
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Both  of  these  restrictions  can  be  summed  up  by  saying  that  the  load 
and  terminal  network  must  remain  electrically  small. 

To  illustrate  the  correlation  between  total  load  current  and 
differential  mode  terminal  current,  WRSMOM  is  used.  Consider  the 
transmission  line  structure  used  for  current  predictions  in  PLOTS 
5-8,  5-9  and  5-10.  These  plots  compare  the  differential  mode  termi¬ 
nal  currents  supplied  by  WRSMOM,  OSMOM,  and  the  transmission  line 
model.  Using  the  same  structure,  the  transmission  line  model  results 
can  be  compared  to  the  total  load  current,  supplied  by  WRSMOM.  For 
this  structure,  WRSMOM  models  the  load  as  one  segment  with  a  match- 
point  at  the  center.  PLOTS  6-1,  6-2  and  6-3  show  the  transmission 
line  model  prediction  for  differential  mode  terminal  current  and  the 
WRSMOM  prediction  of  total  load  current  for  the  ENDFIRE,  SIDEFIRE  and 
BROADSIDE  cases  respectively. 

There  are  several  characteristics  to  note  about  these  plots. 

For  all  three  cases,  the  greatest  deviation  occurs  at  higher  frequen¬ 
cies.  This  is  due  to  the  fact  that  the  load  is  not  behaving  as  a 
lumped  parameter.  (Results  should  be  considerably  better  for  the  wire 
separation  considered  in  Chapter  4,  but  valid  WRSMOM  results  are  not 
available  for  those  cases.)  Unfortunately,  a  great  deal  of  the  devia¬ 
tion  noticed  here  is  not  related  to  the  present  issue,  but  to  the 
pulse  expansion,  point  matching  technique.  As  was  discussed  in  Chap- 
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PLOT  6-1  LINE  DIMENSIONS  (HETERS) 
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PLOT  6“2  LINE  DIMENSIONS  (METERS) 

LENGTH  -  1.0 

EXCITATION:  SIOEFIRE  SEPARATION  «  0.05 


FREQUENCY 


3QD1 1 


ter  3,  this  technique  does  not  allow  adequate  sampling  of  the  E 
field  in  the  termination.  More  zones  are  needed  to  provide  for  more 
matchpoints.  Unfortunately,  adding  more  zones  creates  problems  with 
consistency  of  subsection  length  and  creates  an  overly  large  system 
of  equations.  The  inadequate  representation  of  the  E  field  is  par¬ 
ticularly  important  for  the  ENDFIRE  case,  where  the  method  of  moments 
incorporates  the  incident  field  by  adding  sources  only  to  the  load 
segments . 

Nevertheless,  the  plots  do  illustrate  that  for  this  structure 
and  at  least  for  frequencies  up  through  300  megahertz,  the  common 
mode  currents  do  cancel  in  the  load.  On  this  basis,  the  results  of 
all  three  cases  are  good.  This  holds  true  even  for  the  SIDEFIRE 
case,  where  the  common  mode  current  along  the  line  is  large. 
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CHAPTER  7 


SUMMARY  AND  CONCLUSIONS 

In  any  situation  where  electronic  devices  communicate  with  one 
another  via  transmission  lines,  we  are  faced  with  the  possibility  of 
unwanted  electromagnetic  coupling.  No  matter  what  the  source  may  be, 
the  presence  of  an  incident  electromagnetic  field  can  interfere  with 
this  communication  process.  Given  that  we  must  contend  with  this 
problem,  a  reliable,  accurate  and  efficient  method  is  needed  to  pre¬ 
dict  the  effects  of  this  incident  field.  Transmission  line  theory  is 
an  obvious  choice. 

The  transmission  line  model  offers  several  advantages  when  com¬ 
pared  with  other  methods.  First,  it  probably  provides  the  most  ef¬ 
ficient  solution  technique  known.  This  is  important  when  considering 
several  cases,  such  as  several  different  values  of  loading,  or  solving 
over  a  range  of  frequencies.  Several  thousand  dollars  in  computer 
time  have  been  invested  in  the  work  presented  here;  transmission 
line  model  solutions  made  up  less  than  two  percent  of  that  cost.  A 
second  advantage  to  the  transmission  line  model  is  that  the  mathemati¬ 
cal  equations  involved  are  not  totally  disassociated  from  physical  in¬ 
tuition.  Specific  forms,  allowing  insight  into  the  solution,  can  be 
derived,  as  was  illustrated  several  times  in  Chapter  4.  Transmission 
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line  theory  acquires  these  characteristics  by  making  several  key  as¬ 
sumptions.  These  assumptions  set  forth  certain  requirements  which 
must  be  met  by  the  problem  at  hand.  It  is  the  lack  of  specific  values 
for  these  restrictions  that  prevent  the  establishment  of  limits  of 
validity  for  transmission  line  theory.  It  is  certainly  impossible  to 
validate  the  transmission  line  model  on  a  case  by  case  basis.  The 
results  presented  in  this  work  are  not  intended  as  an  attempt  to  vali¬ 
date  transmission  line  theory  for  all  possible  cases.  Rather,  these 
results  should  be  used  to  indicate  some  of  the  circumstances  where 
and  why  this  theory  may  become  invalid.  At  the  same  time,  it  is  shown 
that  the  transmission  line  model  provides  valid  and  accurate  solutions 
for  several  important  cases  within  its  domain,  when  compared  to  the 
solution  provided  by  the  much  more  rigorous  method  of  moments. 

The  series  of  plots  in  Chapter  4  covers  three  important  loading 
conditions  (low  impedance,  high  impedance  and  matched  impedance)  which 
are  somewhat  representative  of  the  wide  variety  of  possible  loads. 

Each  of  these  loading  conditions  is  examined  for  the  three  most  sig¬ 
nificant  (from  a  transmission  line  theory  standpoint)  cases  of  inci¬ 
dent  field  excitation,  over  a  wide  range  of  frequencies.  The  elec¬ 
trical  separation  between  the  transmission  line  wires  is  small,  never 
becoming  larger  than  1/30  of  a  wavelength  for  all  frequencies  con¬ 
sidered.  In  each  plot,  results  produced  by  the  transmission  line  model 
solution  for  the  terminal  currents  are  compared  with  the  method  of 
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moments  solution  for  the  differential  mode  terminal  currents  pro¬ 
vided  by  OSMOM  (method  of  moments  computer  program  described  in 
Chapter  3) . 

Excellent  correspondence  is  observed  between  the  two  solutions 
for  the  matched  and  high  impedance  load  cases.  The  transmission  line 
solution  remains  accurate  through  at  least  900  megahertz  for  these 
two  loading  conditions.  The  method  of  moments  solution  begins  to 
break  down  at  frequencies  just  below  10  megahertz  for  all  loads. 

This  is  expected  and  due  to  numerical  instability  encountered  during 
matrix  operations.  Correspondence  between  the  two  solutions  did 
not  appear  to  be  significantly  effected  by  the  angle  of  incidence 
or  polarization  of  the  plane  wa\e  used  for  excitation  (at  least  for 
the  three  cases  of  excitation  considered  here) .  The  parameter  which 
seemed  to  have  the  most  effect  in  the  comparison  between  solutions 
was  the  loading,  in  particular  the  low  impedance  loads. 

Although  the  matched  and  high  impedance  transmission  line  model 
solutions  compared  very  well  with  the  method  of  moments,  the  solutions 
for  the  low  impedance  ca.>e  showed  some  slight  discrepancies  at  fre¬ 
quencies  above  300  megahertz.  This  difference  in  solution  is  most 
prominent  at  frequencies  near  those  producing  terminal  current  nulls. 
This  effect  is  similar  to  that  noticed  in  [16]  for  a  transmission  line 
with  short  circuit  terminations.  The  fact  that  this  occurrence  is 
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brought  on  by  low  impedance  loads  might  indicate  that  the  problem  lies 
in  reradiation  by  the  terminal  wires.  The  transmission  line  model 
views  the  load  as  a  lumped  parameter,  thereby  having  no  electromag¬ 
netic  properties.  The  method  of  moments,  however,  views  the  load  as 
a  short  terminal  wire  containing  a  load.  Therefore,  significant 
current  flowing  on  this  terminal  wire  would  cause  it  to  radiate,  as 
an  antenna,  an  effect  which  is  accounted  for  in  the  method  of  moments 
solution,  but  not  in  transmission  line  theory.  Low  impedance  loading 
tends  toward  greater  reradiation. 

To  investigate  reradiation  as  a  possible  cause  for  the  discrepan¬ 
cy,  the  structure  is  slightly  modified  in  Chapter  5.  If  reradiation 
is  the  problem,  lengthening  the  terminal  wire  should  cause  the  cor¬ 
respondence  between  the  two  solutions  to  worsen.  The  wire  separation 
is  increased,  and  indeed  the  correspondence  does  deteriorate  with 
low  impedance  loads. 

Although  reradiation  is  probably  causing  the  greater  part  of  the 
discrepancy,  the  transmission  line  model  begins  to  break  down  for 
other  reasons  for  wide  wire  separation.  As  the  separation  between 
the  transmission  line  wires  is  increased,  the  TEM  assumption,  essen¬ 
tial  to  transmission  line  theory,  begins  to  break  down  and  higher  or¬ 
der  modes  become  significant  as  they  begin  to  propagate.  To  see  the 
effects  of  this  assumption  losing  its  validity,  reradiation  must  be 


-139- 


removed.  To  accomplish  this,  the  low  impedance  loads  are  replaced  by 
open  circuits.  Since  this  will  obviously  cause  the  terminal  currents 
to  become  zero,  the  two  methods  of  solution  can  be  compared  by  looking 
at  various  points  along  the  line.  The  transmission  line  model,  via 
the  chain  parameter  matrix,  can  generate  the  differential  mode  cur¬ 
rent  at  any  point  along  the  line.  The  method  of  moments,  by  the  na¬ 
ture  of  the  technique,  produces  current  values  at  various  points  on 
the  structure.  This  method  of  moments  solution  for  total  current  can 
then  be  decomposed  into  common  and  differential  mode  currents,  and 
the  differential  mode  current  can  be  compared  to  the  transmission  line 
model  solution  for  these  points  along  the  line.  Three  separations  are 
considered.  Each  time  the  separation  is  increased,  the  transmission 
line  solution  deviates  further  from  the  method  of  moments  solution,  as 
expected.  For  the  widest  separation,  the  ratio  of  line  length  to  sepa¬ 
ration  is  10:1 ,  causing  a  questionable  transmission  line  model  solution. 

Transmission  line  theory,  by  the  TEM  assumption,  assumes  that 
common  mode  currents  are  not  present.  That  is,  it  solves  only  for 
differential  mode  currents.  This  raises  a  significant  question.  Of 
what  value  are  the  transmission  line  results  if  they  only  predict  the 
differential  mode  component  of  the  total  current?  Remembering  that 
transmission  line  theory  views  the  load  as  a  lumped  parameter,  the 
common  mode  currents  cancel  in  the  load.  Therefore,  only  the  dif¬ 
ferential  current  flows  through  the  load.  This  remains  true  as  long 
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as  the  assumption  that  the  load  behaves  as  a  lumped  parameter  remains 
valid;  the  load  must  remain  electrically  small.  High  frequency  and 
wide  separation  will  cause  the  use  of  the  transmission  line  model 
solution  as  an  accurate  prediction  for  total  load  current  to  be 
questionable.  In  Chapter  6,  the  transmission  line  model  solution  is 
compared  to  the  total  load  current  prediction  provided  by  WRSMOM 
(method  of  moments  computer  program  discussed  in  Chapter  3) .  The 
results  show  that,  for  this  separation  the  transmission  line  model 
solution  provides  an  accurate  prediction  of  total  load  current  up 
through  300  megahertz,  or  when  the  separation  is  1/20  of  a  wavelength. 

At  this  point,  a  few  comments  should  be  made  concerning  the 
method  of  moments  computer  program  used.  0SM0M  [13],  [14],  which 
uses  Galerkin's  method  and  piecewise  sinusoidal  expansion  functions, 
was  very  effective  in  producing  baseline  data  for  comparison.  It 
seems  to  handle  the  comers  and  terminations  very  well  and  appeared 
to  give  very  reliable  results  for  all  cases  considered.  As  was  dis¬ 
cussed  in  Chapter  5,  some  problems  were  encountered  with  the  pulse 
expansion  and  point  matching  technique  used  in  WRSMOM  [12] .  Although 
there  appear  to  be  no  programming  problems,  the  technique  itself 
does  not  seem  to  be  very  well  suited  for  transmission  line  analysis. 

It  often  becomes  difficult  to  model  a  transmission  line  structure 
and  maintain  the  consistency  in  subsection  length  necessary  for  a 
valid  solution.  Also,  in  many  cases  the  E  field  cannot  be  adequately 
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sampled  in  the  termination  to  provide  a  valid  solution,  using  point 
matching.  This  is  discussed  in  further  detail  in  Chapters  3  and  5. 

In  general,  the  results  of  this  investigation  show  that  the 
transmission  line  model  can  be  successfully  applied  to  a  certain 
class  of  problems.  The  efficiency,  intuitive  nature  of  the  equations 
and  accuracy  in  the  lower  frequency  range  make  it  particularly  invit¬ 
ing  as  a  solution  technique.  The  user  must  be  conscientious  in  the 
observance  of  restrictions  placed  on  the  method  by  the  assumptions  of 
transmission  line  theory.  The  transmission  line  model  is  obviously 
not  designed  to  solve  all  types  of  coupling  problems.  However,  it 
does  appear  to  serve  a  class  of  problems  commonly  encountered  in 
present  wiring  techniques.  The  next  step  in  the  process  of  estab¬ 
lishing  validity  to  the  technique  is  to  further  define  the  restric¬ 
tions  on  the  electrical  dimensions  of  the  structure  (restrictions 
required  for  a  valid  transmission  line  model  solution),  e.g.:  wire 
separation  defined  in  terms  of  wavelength.  This  step  will  probably 
be  difficult  because,  as  shown  earlier,  several  factors  intervene, 
e.g.:  low  impedance  loading.  Used  correctly,  the  transmission  line 
model  can  be  a  very  effective  and  useful  solution  technique. 
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The  use  of  an  incident  electromagnetic  plane  wave  to  excite  a 
two  wire  transmission  line  is  described  in  this  appendix.  The  plane 
wave  is  assumed  to  be  traveling  in  free  space  away  from  the  origin 
of  the  coordinate  system.  Four  parameters  are  necessary  to  describe 
the  wave.  Two  angles  are  used  to  describe  the  direction  of  propa¬ 
gation.  (See  FIGURE  1-1.)  The  angles,  9p  and  <j>p,  are  very  similar 
to  the  angles,  9  and  4> ,  encountered  in  a  spherical  coordinate  system. 
9p  is  the  angle  between  the  y  axis  and  the  propagation  path  of  the  in¬ 
coming  wave.  (j>p  is  measured  from  the  z  axis  to  the  projection  of  the 
propagation  vector  in  the  x-z  plane.  A  third  angle,  9£,  is  needed  to 
describe  the  polarization  of  the  plane  wave  (the  direction  of  the  E 
field).  Assume  a  plane,  containing  the  origin,  and  perpendicular  to 
the  propagation  vector.  Because  of  the  TEM  structure  of  the  wave, 
the  E  field  vector  will  lie  in  this  plane.  9p  is  measured  clockwise 
(looking  in  the  direction  of  propagation)  from  the  projection  of  the 
y  axis  (onto  the  plane)  to  the  E  field  vector.  The  last  parameter 
needed  to  describe  the  wave  is  the  magnitude  of  the  E  field.  The 
H  field  magnitude  can  be  related  to  the  E  field  magnitude  by  the  in¬ 
trinsic  impedance  of  free  space,  n- 
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where 


n  =  ^7^  (1-2) 

and  yQ  is  the  permeability  of  free  space  and  eQ  is  the  permittivity 
of  free  space.  The  orientation  of  the  H  field  can  be  obtained  by 
equation  (2-52). 

The  plane  wave  is  propagating  with  a  propagation  constant  of  k. 

For  free  space  k=B,  described  in  Chapter  2.  Having  the  plane  wave 
so  defined  still  leaves  us  with  the  problem  of  obtaining  the  components 
of  k  and  Em  (E  field  magnitude)  in  x,  y,  and  z  for  use  in  the  trans¬ 
mission  line  equations.  This  becomes,  for  the  most  part,  a  coordinate 
transformation.  The  three  components  of  k  can  be  obtained  by  a  trans¬ 
formation  from  spherical  to  rectangular  coordinates,  as  described  in 
[5].  Given  the  magnitude  of  k,  the  following  equations  yield  values 
for  the  rectangular  components  of  k. 

kx  =  k  sin9p  sin<|>p  (1-3) 

ky.  =  k  cos 0p  (1-4) 

kz  =  k  sin6p  cos$p  (1-5) 
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Obtaining  the  rectangular  components  of  the  E  field  is  slightly 

more  difficult.  Given  the  value  of  E  ,  the  total  E  field  magnitude, 

m 

the  values  of  the  rectangular  components  are  as  follows: 

E  =  -  E  cos9_  cos6n  sinc>_  -  E  sin8_  cos$n  (1-6) 
xm  m  E  P  P  m  E  P 

E  =  E  cos0_  sin0n  (1-7) 

ym  m  E  P 

E  =  -  E  cos6_  cose,,  +  E  sin0_  sin<t>_  (1-8) 

ZIH  C  r  C  r 
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APPENDIX  II 


In  Appendix  III  a  listing  is  given  for  the  computer  program, 
written  in  FORTRAN,  which  was  used  to  predict  the  terminal  currents 
of  the  two  wire  structure,  by  transmission  line  theory  techniques. 

The  format  of  the  input  data  is  also  given.  This  appendix  describes 
that  program  and  the  required  input  data. 

The  first  task  of  the  program  is  to  define  the  wire  structure 
using  the  parameters  discussed  in  Chapter  2.  The  data  inputs  necessary 
are  wire  radius,  wire  separation,  wire  length,  the  load  at  x=0,  and 
the  load  at  x=L.  The  program  will  calculate  the  transmission  line 
solution  for  arbitrary  loads;  that  is,  there  are  no  constraints  on 
the  loads.  The  loads  at  x=0  and  x=L  are  independent  (they  can  be 
different  from  each  other)  and  can  be  complex.  The  next  inputs  de¬ 
scribe  the  uniform  plane  wave  used  as  the  incident  field  excitation. 

The  four  parameters  necessary  to  define  the  wave  are  Em,  0  ,  9p,  and 
4>p  as  described  in  Appendix  I. 

All  calculations  are  done  using  double  precision  to  assure  ac¬ 
curacy.  A  special  effort  was  made  to  do  all  calculations  using  a 
numerically  stable  form  of  the  equations  involved.  For  example,  the 
calculation  of 
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(II-l) 


e 


First,  calculations  that  do  not  involve  frequency  are  performed. 
The  rectangular  coordinate  components  for  the  E  field  are  obtained  as 
described  in  Appendix  I.  Next,  the  per-unit-length  capacitance  is 
calculated  using  the  wire  radius  and  separation,  as  given  in  [6] . 


Once  the  value  for  c  has  been  obtained,  Rc>  the  characteristic  resis¬ 


tance  of  the  line  is  calculated  as  given  by  equation  (2-50) . 


(H-4) 


The  next  section  performs  operations  involving  frequency  for  each 
of  the  frequencies  input  as  data.  The  calculations  start  with  ob¬ 
taining  8  and  rectangular  components  of  k.  The  calculation  of  1(0) 
and  I  CL)  is  done  using  slightly  modified  versions  of  the  forms  pre- 
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APPENDIX  III 


The  following  three  pages  contain  a  listing  of  the  computer  pro¬ 
gram  which  was  used  to  generate  all  of  the  transmission  line  model 
results  presented  in  Chapters  4,  5,  and  6.  A  description  of  this  pro¬ 
gram  is  presented  in  Appendix  II.  The  format  of  the  input  data  is 
shown  on  the  page  immediately  following  the  program  listing.  Multiple 
runs,  for  different  frequencies,  can  easily  be  executed  by  simply 
adding  the  desired  frequencies  to  the  end  of  the  data  string,  using 
the  format  shown.  The  program  will  produce  a  solution  for  each  fre¬ 
quency  and  stop  execution  when  it  runs  out  of  data  (when  it  runs  out 
of  frequency  data  entries) . 

Output  produced  by  this  program  was  compared  to  output  pro¬ 
duced  by  the  program  WIRE  (17)  for  several  representative  test  cases. 
The  results  from  both  programs  were  the  same. 
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IMPLICIT  REAL*8 (A ,  B ,  D  ->H , 0 -Z)  , C0MPLEX*1 6(C)  00000010 
REAL* 8  LEN ,L  00000020 
CJ«DCMPLX(0.D0,1 .DO)  00000030 
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APPENDIX  IV 


also 


wvic  =  ojv(*€c)^  =  wv(i)  =  ^  =  k 


Making  these  substitutions. 


!'( L)  =  jkE  d  f  [cos (k(L-x))  +  j  sin(k(L-t))]  e'^dx 
s  m  Jo 

.  jkEd  fL  ejk(1-T)  e'ikTd, 

=  jkE  d  f  e  ^kTdx 

ra  Jo 


=  jkE  d  e 
J  m 


jkL 


-  j  2kx 


"=PF 


m,  jkL  ,  -j2kL 
*  -  yd  eJ  (e  J  -  1) 


/e'jkL  .  e*kL\ 

■  *  VC— r1-- 


V'(L)  =  jE  d  sin(kL) 
s  m  1 


(IV-4) 


Substituting  equation  (IV- 13  and  (IV-2)  into  equation  (2-43) 
gives  equation  (IV-5).  The  steps  of  evaluation  leading  to  a  simpli¬ 
fied  form  for  I'(L)  follow. 


IsCL)  *  /  I-  jvc  sin (k(L-i)jV  (T) 
0  s 


+  cOS(k(L-T))I oCT)JdT 


/L- 

J  n 


=  E 


Note  that 


jvc  sin(k(L-T))  jaii  J“d  e'jkT 
o  n 


-  cosCk(L-r))  jaicE^d  e’jkx]dT 

f  L 

j  I  vcaiij 

1  J0  n  slnCk(L-T))  -  juc  cosCk(L-r))  e'^Tdj 


uy  _  u  __i_  _  i 

n '  =vWT^ i  =  1 


(IV-5) 


r L 

*  -  j-cE^d  (cos(k(L-t))  .  j  sin(k(L-T))]e-ikTl 

f  L  ejk<L-' 

J  0 

-  juicE  d  / 

m  Jo 


*  -  ja.cE  d  /  1  e^k(L-T)  e'jkTdx 


'  eJ'kL  e-j^r 


,-j  2kT 


-i  L 


5JE 


i  d  eJ 
m 

JkL  , 
e'  (e 

-j  2kL 

e'JkL  -  ejkL\ 

1) 


aicE  d 

-  J  — 1~  sin(kl) 


-157- 


VMxsm*..  ■-  .  -^%1-nti  iiiiM^^i- 


Note  that 


uc 
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■  VC  ■ 
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Ed 

I/1)  ■  -  j  f  sin(kL) 


(IV-6) 
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APPENDIX  V 


The  purpose  of  this  appendix  is  to  show  that  for  ENDFIRE  and 
BROADSIDE  excitation,  with  open  circuit  loading,  the  transmission  line 
solution  of  Chapter  2  gives  zero  currents  for  all  values  of  x  along 
the  line. 

Equation  (5-6)  has  been  specialized  to  open  circuit  loading  but 
is  still  valid  for  all  three  cases  of  excitation.  It  is  given  here, 
as  equation  (V-l). 


(v-i> 

First  consider  the  ENDFIRE  case.  The  value  of  I'(L)  is  given  by 
equation  (4-9). 

E  d 

Is(L)  =  ’  jir  sin(eL)  (v-2) 


Variable  substitution  yields 


Ed 

rs(x)  =  '  jf  sin  (ex)  (V-3) 


Substituting  (V-2)  and  (V-3)  into  (V-l),  we  obtain 


i** 
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(V-4) 


.  -..Ed  Ed 

T,  ..  .sm(Bx)  m  .  ....  .  m  .  . 

I(x)  *  TT  sin^L>  -  J TT  sln ^ 

c  c 

Simplification  yields  the  following. 

Ed  Ed 

1 00  =  j*gr-  sin(Bx)  -  sin(Bx)  =  0  (V-5) 

c  c 

the  BROADSIDE  case  is  very  similar.  The  value  of  I'(L)  fspecialized 
to  the  BROADSIDE  case) ,  as  given  by  equation  (4-38) : 

E  d 

Is(L)  “  '  jir  sin(6L) 
c 

This  relation  is  identical  to  equation  (V-2).  Therefore  the  result 
will  be  exactly  the  same  as  that  given  by  equation  (V-5). 
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